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Summary. — In the past few years the field of hadron spectroscopy has seen 
renewed interest due to the pubblication, initially mostly from B-Factories, of evi- 
dences of states that do not match regular spectroscopy, but are rather candidates 
for bound states with additional quarks or gluons. A huge effort in understanding 
the nature of this new states and in building a new spectroscopy is ongoing. This 
report reviews the experimental and theoretical state of the art on heavy quarko- 
nium exotic spectroscopy, with particular attention on the steps towards a global 
picture. 

PACS 14 . 40 . Rt - Exotic mesons. 
PACS 14.40.Pq - Heavy quarkonia. 
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1. Introduction 

The discovery of the J/ip in 1974 has been of pivotal importance for the development 
of the Standard Model of particles. On the one hand its existence confirmed the GIM [1] 
prediction of a fourth quark, the charm. With the addition of the charm to the u, d, s 
quarks, the problem of flavor changing neutral currents in the Cabibbo theory was re- 
moved and a 2 x 2 mixing matrix of quarks was found for the first time. On the other 
hand the impressive narrowness of the J/ip (about 90 keV) was realized to be the foot- 
print of the asymptotic freedom of quantum chromodynamics. Thus the J/ip discovery 
opened a new route to the theory of weak interactions and pointed to the most striking 
feature of the theory of strong interactions, namely the quark freedom on very short 
timescales. 

The bound state of two charm quarks, usually named charmonium, is the strong 
interaction analog of an atomic system. The reduced mass system of the two charm 
quarks is subject to a central potential of the form V(r) = a\/r + a 2 r, where the first 
term is the Coulomb-like potential between color charges and the linear term is the 
phcnomcnological implementation of the confining force between quarks. 

In the rcnormalization of quantum chromodynamics, a scale, A QCD , arises by 'dimen- 
sional transmutation' of the ultraviolet cutoff. In the pure gauge theory, A QCD enters in 
the definition of the running coupling constant a s (/x) = 27r/[ll ln(^,/A QCD )] showing that 
perturbation theory is valid only at distances 1/fi smaller than 1 / A QCD . A priori A QCD is 
an arbitrary scale which can take any value. Taken over the full set of available data (from 
deep inelastic scattering, T decay, e + e~ — > jets etc) one finds 100 < A QCD < 400 MeV. 
The fact that it turns out to be A QCD > > m u , is at the heart of the reason why light 
quarks are confined in hadrons of the 1 fm size. 

On dimensional grounds we can state that a restoring force between quarks could be 
written as F q q ~ A^ cd , thus resulting from a potential V(r) <~ A^ co r. Fitting a 2 on 
experimental data one gets a 2 = I /a 2 with a = 2.3 GcV -1 . Observing that a 2 ~ A^ CD 
one finds indeed A QCD ~ 400 MeV. This in turn means that the effects of confinement are 
apparent at a distance 1/A QCD = 0.5 fm. But the J/ip is three times heavier, and thus 
smaller than the proton. Thus one can imagine that the quarks in the J/ip are so close, 
i.e. far from the confinement distance 1/A QCD , to be considered as almost free. The 
mass of the charm m c is therefore just half of the mass of the J/ip, m c ~ 1530 MeV, and 
non-relativistic methods can be used to describe charmonia. In this respect the discovery 
of the J/ip can be considered as the first evidence of the existence of quarks. 

Since the discovery of the J/ip, the charmonium system has been thoroughly stud- 
ied [2]. Very precise determinations of spectra have been provided including hyperfine 
effects. The decay mechanisms have been investigated in greater detail and no one ex- 
pected surprises from charmonia. Particles like the rj c and predicted and searched 
since many years, have eventually been detected, with no striking unexpected features. 

It was therefore particularly impressive to discover in 2003 a charmonium-like system, 
the V(3872), decaying into J/ipmr, thus behaving as a charmonium excitation, but 
surprisingly with almost equal decay rate into J/ip p and J/ipoj. Not to speak of the 
fact that the V(3872) was soon realized not to have the very well known radiative decay 
pattern of the closest in mass standard charmonium state: the radial excitation 2 1 P 2 . 

After the discovery of the V(3872) more and more similar narrow resonances have 
been discovered and confirmed at electron-positron and proton-antiproton colliders. Namely, 
since the discovery of the V(3872), about twenty new unexpected charmonium-like par- 
ticles have been found. Most of them have definitely clear clashes with standard char- 
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monium interpretations. 

These facts opened a debate about the nature of the newly discovered particles. The 
mainstream thought has been that of identifying most of these resonances as molecules of 
open charm mesons. In particular the A (3872) happens to have a mass very close to the 
DD* threshold. The binding energy left for the X is consistent with £ <~ 0.25±0.40 MeV, 
thus making this state very large in size: order of ten times bigger than the typical range 
of strong interactions. Relying on intuition, this feature, per se, would discourage to 
pursue the molecular picture of the X: how is that possible that such a loosely bound 
state can rearrange its quarks to produce a J/tpp final state? It must occur that far 
apart charm quarks must coalesce in a very compact object as we described the J/tp to 
be. How is that possible that the X has such a large prompt production cross section 
at CDF as 30 -j- 70 nb? In a high energy collision environment one can think that it is 
rather difficult to form a molecule with practically zero binding energy. Why is the decay 
width of the X in DDn not exactly equal to the known decay width of the D* meson 
component of the state? In particular what would be accelerating the decay rate of the 
D* , as data seem to suggest, inside an almost unbound state? 

These questions apply to other near-to-threshold hypothetical molecules and have 
induced to think to alternative explanations for the X and its relatives. A rather natural, 
and more fundamental possibility, would be that of thinking to other forms of aggregation 
of quarks in hadrons, like diquarks. A diquark is a colored quark-quark state which could 
neutralize its color, binding with an antidiquark. The resulting meson is a particular 
realization of a tetraquark. Indeed as for the color, the diquark is like an antiquark and 
the antidiquark is just like a quark. It follows that the resulting tetraquarks are kind of 
standard qq mesons but with the notable difference that at the end of the electric color 
string there are diquarks instead of quarks. 

What would happen then when stretching the color string at a distance 1 //x larger than 
1/A QCD ? The standard lore for conventional mesons is that the string, which provides 
the restoring force F qq -, would break at some point. By then it stored enough energy 
to excite the vacuum to produce a pair of light quarks at the breaking point. These 
in turn saturate the color of the string endpoint quarks, turning the decaying meson 
into two lighter mesons. This pictorial way to describe meson decays has a quantitative 
implementation in all the hadronization algorithms used in Monte Carlo simulations of 
collider physics events. 

Following the same line of reasoning, a tetraquark is expected to decay into two 
baryons upon color string breaking. Indeed diquark-quark configurations are possible in 
standard baryons as shown by a very simple interpretation of the ratio of deep inelastic 
scattering (DIS) neutron and proton structure functions in the x — > 1 limit. Indeed it 
is known that experimentally the ratio (a^/i^ (x) — > 1/4 as the fraction x of the 
momentum carried by the parton involved in DIS tends to 1. It results that jt^ 11110 ' 0011 ) — 
^l2q e qfq( x )i where f q [$ the parton distribution function for the q quark species. As 
the quark participating to DIS gets closer to carry the entire nucleon momentum, the 
remaining quarks are frozen in their lowest energy state. If diquark bound states exist, 
these will further lower the energy. But, as shown by lattice simulations, diquarks of 
the form [uu] or [dd] cannot be formed. Thus the lowest energy configurations for the 
spectator quarks in DIS are reached when we have a [ud] diquark in the neutron with a 
d quark involved in DIS and again a [ud] diquark in the proton with a u quark in DIS. 
Since e d = -1/3 and e u = 2/3, the experimental fact F^ n) (x) / F^ p) (x) -> 1/4 (a; 1/4) 
is understood as due to the ratio of charges e^/e^ = 1/4. 
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The X(3872) on the other hand has not enough mass to decay into two charmed 
baryons because of phase space. This forces the diquark-antidiquark system to rearrange 
itself into a J/ip p or J/ip lo configuration. Contrary to the picture given in molecular 
models, such a rearrangement happens inside the boundaries of 1/ A QCD . But some of the 
newly discovered hadrons have enough mass to decay into two baryons. The F(4660), for 
example, appears to decay prominently into A C A C baryons, as expected for a tetraquark. 

The tetraquark model has more challenging predictions though. Charged states are 
expected, such as [cm] [erf] and even doubly charged states as [cu][Js]. The model does not 
predict light doubly charged states though (since diquarks like [uu] would be involved). 
In particular, in one of the first papers on the tetraquark model [3], it was predicted 
that narrow resonances decaying into J/ipn + should be observed. These kind of particles 
have indeed been observed by Belle about three years later. 

The Z+(4430), the first of a series of three newly discovered charged particles, decays 
in a charmonium state and a charged pion. Even if Belle observes it at more than 6a 
significance, BaBar has only 3.5<r and CDF has not yet confirmed it. Moreover the 
tetraquark model predicts charged partners of the (3872) with a very similar mass. 
The latter have not yet been observed. The doubly charged states, as for the time being, 
have not been searched. 

In order to explain in the tetraquark model the peculiar decay pattern of the X(3872) 
which dissociates with equal rate into J/ ipp and J/ipui one needs to consider two neutral 
states with a difference in mass of few MeV. One can call them X u = [cu][cu] and 
Xd = [cd] [cd] . Once these states are mixed by uu -R- dd annihilations we can end up with 
isospin pure mixtures of the kind (X u ± Xd)/V2. But if we suppose that annihilation 
contributions are suppressed by the smallness of a s at the charm mass, we could align 
the eigenvectors of the mass matrix in the quark basis rather than in the isospin basis. 
In that extreme case we have two unmixed states X u — [cu] [cu] and Xd = [cd] [cd] : each 
of them can contain both 1 = and 1=1. The difference in mass between the two must 
be Am <~ nid — m u . As of today it seems that the 3.5a evidence of both a X(3872) and 
a X(3875) resonances have faded out by Belle and BaBar data analyses. (At some stage 
we thought that the former was decaying preferably in DD* and the latter in J/ipp. The 
decay pattern of these two neutral states was explained in [4]). 

The main drawback of the tetraquark picture, as it will be discussed in this review, 
is the proliferation of expected states. We do not have any clue of selection rules which 
could limit the production of tetraquark particles at the fragmentation level. Therefore, 
as it will be shown in the following, we will merely enumerate all the possible bound 
states trying an estimate of their mass values. Unfortunately we cannot predict the 
fragmentation probabilities of tetraquark states, so it is also difficult to estimate their 
production rates. 

The mass spectra are computed relying on the non-relativistic constituent quark 
model where the quark masses (constituent masses) incorporate the bulk of the non 
perturbative interactions and basically are fitted from baryon and meson spectra. Cor- 
rections to the simple algebraic sums of those masses come from chromomagnctic (spin- 
spin) interactions among constituents. Since we will be dealing with diquark-antidiquark 
mesons, spin-spin interactions will be of two kinds: those occurring inside single diquarks 
and those between the two diquarks. As the relative orbital angular momentum of the 
tetraquark grows (L > 0), point-like spin-spin interactions will be switched off. Actually, 
in this review, we propose to treat the orbital excitations and the spin-orbit corrections 
in the framework of a relativistic string model. The resulting model has been table tested 
on standard charmonia and has a rather good behavior at determining mass spectra. 
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One should be aware that the 'realm' of the constituent quark model is actually the 
one of very small A QCD , actually A QCD << m u ,md- In this hypothetical world the pion 
would weight 600 MeV and there would be a clear clash between the Nambu-Goldstone 
interpretation of light pseudoscalar mesons and their constituent quark nature. Yet, it 
is surprising to observe that the constituent quark model is a powerful tool in hadron 
spectroscopy. An explanation of this from first principles likely requires new ideas in 
non-perturbative chromodynamics. 

The scope of this short review is to give a summary of the tetraquark model expec- 
tations and of the methods used to compute masses and decays in this framework. We 
will briefly discuss the molecular interpretations of the X, Y, Z mesons, and give a com- 
prehensive overview of the experimental situation. We will clearly stress the drawbacks 
of the various models. In particular, in our discussion we will show that the (3872) has 
serious difficulties both with the tetraquark and the molecule interpretation. 

We aim at providing a thorough discussion of what has been made and what is still 
waiting to be done on the experimental side in order to help the understanding of the 
new hadrons. 

2. Theoretical Models 

Quark, anti-quarks and gluons do not necessarily bind only in the qq pairs or qqq 
triplets. The quark model from its beginning [5] contemplate the existence of other 
aggregations. The states recently discovered suggest in particular the possibility to be 
observing states with two or three quarks and two or three anti-quarks or a quark, an 
anti-quark and valence gluons. While the latter states (called hybrid) appear in only one 
configurations, there are two possibilities to form bound states out of two quarks and 
two anti-quarks: 

• binding the two quarks in a colored configuration called diquark [qq] a , ot anti- 
diquark [qq] a . Color charge is neutralized by the interaction diquark- antidiquark: 
this configuration is what will be called tetraquark; 

• binding each quark to an anti-quark [q a q a ] and allowing interaction between the 
two color neutral pairs [g a <f*] [g.ag' 3 ]. This configuration is very close to the one 
with two interacting mesons and it is usually called molecular bound state. 

This section details the differences between the approaches in terms of predictions 
and interpretations of the currently observed states. Besides the exotic interpretations 
the possibility is still open of at least a fraction of them being ordinary charmonium 
with an incorrect J assignment or for which the potential model predictions do not 
hold. It is to be considered indeed that, as shown in Fig. 6, while below the open charm 
threshold all the narrow states have been observed and match predictions, above it almost 
all conventional charmonium states are missing and the few candidate ones have masses 
significantly different from predictions. Another lurking possibility is that the observed 
states are not strong resonances but the effect of opening thresholds. For a discussion of 
this aspect we remit to other reviews as for instance Refs. [6, 7]. 

Finally we briefly review the picture, recently proposed, in which some of the observed 
states are described as a standard charmonium stuck inside a light hadron, the so-called 
hadrocharmonium. Some hints about pentaquark and hexaquark structures are given in 
the end of this section. 
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2T. Tetraquarks. - Tetraquarks are bound states of a diquark and an anti-diquark, 
combinations of two quarks and two anti-quarks respectively. The combination of two 3 C 
color states can yield a 3 C or a 6 C . Which of these state are bound can be estimated in 
the one gluon exchange model[8],[9]. As in electromagnetism, the discriminant between 
attractive and repulsive forces is the sign of the product of the charges: in the case 
of the color interaction, the sign of the product of SU(3) C generators T a for the A(B) 
representation is computed as 

(1) Joc^TJTg = \J2 [T a2 -Tf -T a A o=l,---,8; A, B = 3 or 3 

a a \ / 

where 

(2) T a = T% + T% 

is the generator on the tensor product A ® B space. Since T° 2 = ^ a T a T a commutes 
with the SU(3) generators it is a number on each irreducible representation (Casimir). 
Using the Shur lemma, the Casimir operator is proportional to the identity I: 

(3) ^T a2 = C(D% 

a 

where C is a constant depending on D, the dimension of the tensor product space gen- 
erated by T a — T% + Tg . The discriminator then becomes 

(4) / oc \{c{A ® B) - C{A) - C{B)^ 
Imposing the normalization of the Casimir operator one obtains 

where 

(6) k A ®B = k A + k B 

kA®B = dim{A)kB + dim(B)kA 

Tab. I shows the estimated kn coefficients as evaluated in the one gluon exchange model 
analysis. 

These results show that only the 3c state can be bound and therefore this is the only 
one considered hereafter. It is to be noted that the same reasoning, applied to a quark 
anti-quark pair, shows that only the color singlet is bound, as known. 

The other relevant quantum number of the diquarks is the spin, which can be S = 
or S = 1. Lattice simulations [10] show that scalar diquarks made of light quarks are 
more stable than vector ones: for this reason they are usually called "good" and "bad" 
diquarks respectively. 

Predictions of the tetraquark model are obtained with three different approaches: 
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Table I.: One gluon excange model results. 



D = A® B 


1 


3 


6 


8 


k D 





1/3 


5/2 


3 


C(D) 





4/3 


10/3 


3 


I=UC(D)-C(A)-C(B)) 


-4/3 


-2/3 


1/3 


1/6 



• QCD sum rules: calculations are based on the correlation functions of two hadronic 
currents with the quantum numbers of the hadron under investigation [11, 12]. At 
the quark level the two-point correlation function is expanded as a series of local 
operators (generalized Wilson OPE): 



(7) i[(q 2 ) = i / tfxe*^0\T\j{x)3\G)W) = £c n (Q 2 )6„, 

J n 

where the coefficients C n (Q 2 )(Q 2 = —q 2 ) include by construction only the short- 
distance domain and can therefore be evaluated perturbatively, while the set {O n } 
include all local gauge invariant operators expressible in terms of gluon and quark 
fields. The fundamental hypothesis of the sum rules approach is that the current- 
current correlation function in Eq. (7)can be matched to the dispersion relation 

(8) U^) = - fds 2 PiS \. +■■■ 
11 J q z - s + is 

To extract physical quantities, the spectral function p(s) is studied by parame- 
terizing it as a single sharp function with a pole representing the lowest physical 
resonance and a smooth continuum representing higher mass states. A modern 
introduction to the method of QCD sum rules can be found in [13]; 

• Relativistic quasipotential model [14, 15]: the interaction between two quarks in a 
diquark and between diquarks and antidiquarks in a tetraquark are described by 
the diquark wave function ^ of the bound quark-quark state and by the tetraquark 
wave function of the bound diquark-antidiquark state respectively. They must 
satisfy the Schrodinger equation: 

(9) - *„(p) - J |^,Hp.<.; «)*„«■), 

where the relativistic reduced mass is: 



(10) 



Ei + E 2 
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and Ei, E 2 are given by: 



_ M 2 - mj + m\ _ M 2 - m 2 + mj 

1 J 1_ 2M ' 2 2M 

Here M = E\ + E 2 is the bound state mass (diquark or tetraquark), m\$ are 
the masses of the quarks which form the diquark or of the diquark and antidi- 
quark which form the tetraquark; p is their relative momentum. In the center of 
mass system the relative momentum squared on mass shell is indicated as b 2 (M). 
The kernel V^tCp, q; M) is the potential operator of the quark-quark or diquark- 
antidiquark interaction and depends on the quark composition of the object under 
investigation. 

The potential given by Eq. (9) is then solved numerically for the diquark wave 
function VE^, the diquark masses and form factors. 

Next, the masses of the tetraquarks are calculated introducing a kernel for the 
diquark-antidiquark bound state that involves the diquark parameters found in the 
first step. 

• A generalization of the constituent quark model: hadron masses are described by 
an effective Hamiltonian that takes as input the constituent quark masses and 
the spin-spin couplings between quarks. By extending this approach to diquark- 
antidiquark bound states it is possible to predict tetraquark mass spectra. This 
model will be discussed in detail in the next section. 

2' 1.1. Spectra. We estimate the tetraquark mass spectra in the framework of the 
non-relativistic constituent quark model combined with a hadron string model. In the 
following section the two models arc defined and the mathematical details for both arc 
given. The numerical inputs required are discussed and derived and the results for the 
mass spectra determination are summarized for different quantum numbers. 

In order to define the energy of a tetraquark state we have to consider all the possible 
interactions between the quarks. In the ground state the two diquarks interact only 
by spin couplings because the angular momentum is zero (L — 0). An effective non- 
relativistic Hamiltonian can be written including spin-spin interactions within a diquark 
and between quarks in different diquarks. 

If the angular momentum is different from zero (L ^ 0) the pointlike spin interactions 
are suppressed as the average distance between diquarks grows. Tetraquark in orbitally 
excited states are studied in a model of color flux tubes first proposed by Selem and 
Wilczek [16] in order to emphasize the importance of diquark correlations in hadronic 
physics. We will generalize the model to describe a four quark state as a spinning string 
with a given angular velocity to connecting two massive diquarks mi, m 2 at distances r\ 
and r 2 away from the rotation axis. 

The non-relativistic approach for L = 

The most general Hamiltonian describing the L = tetraquark state of flavor com- 
position [<7i<j2][9394] can be written as: 

(12) H = mm+mq2 + H ( f ) +H (f 
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Table II.: Numerical values of quark-antiquark spin couplings. 







sq 


ss 


cq 


cs 


cc 


bq 


bs 


be 


bb 


(MeV) 
n qq (MeV) 


315 
103 


195 
64 


121 

22 


70 

25 


72 
72 


59 
6 


23 
8 


23 


20 


36 



where m qi are the diquark masses and: 

(13) H^ q ^ = 2n qi (S qi ■ Sq 2 ) + 2n q2 (Sq 3 ■ S q4 ) 

(14) = 2n qiq3 (S qi ■ Sg 3 ) + 2n qiqi (S qi ■ S q4 ) + 2n q2q3 (S q2 • S q3 ) + 2n q2qi (S q2 • S qi ). 

the first term refers to spin-spin interactions between quarks bound to give a diquark 
and the second term refers to interactions between quarks in different diquarks. The 
coefficients n qiqj depend on the flavor of the constituents qij and on the particular color 
state of the pair. We remind that the coefficients n qiq , include the diquark masses 
dependence and the point-like behaviour of the spin-spin interactions: 

(15) K mqj S qi ■ S qj ~ ~^-S qi ■ S qj 6 3 (fl rj) 

It follows that the coefficients K qiqj have the dimention of energy. The constituent quark 
masses and the K qqi couplings for color singlet combinations are determined from the 
scalar and vector light mesons. They can then be translated into the tetraquark couplings 
using 

(!6) K qiq3 {[qm] [9394]) = ^{K qi q 3 )l 

as detailed in Ref . [3] . The n qq couplings are instead determined from the masses of the 
qqq baryons ground (J = 1/2) and excited (J = 3/2) states. 

The numerical values for the free parameters useful for the determination of tetraquark 
mass spectra are summarized in Tab. II-III. 

To obtain the diquark masses one state must be assumed as tetraquark and the rest 
of the spectrum can be derived accordingly. In this review will start assuming that the 
X(3872) is a {cq]s=i[cq]s=o tetraquark (see Sec. 4T). By diagonalizing the Hamiltonian 
in Eq. (12) and, using the spin couplings values derived above, we obtain the mr cg i 
diquark mass value. In order to reduce the experimental information needed we estimate 
the remaining diquark masses by substituting the costituent quark forming the diquark. 



Table III.: Constituent quark masses as determined from meson. 





q 


s 


c 


b 


m, (MeV) 


305 


490 


1534 


4720 
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Table IV.: 


Diquark masses. 








[cq] 


[cs] 


M 


[bs] 


m qiV (MeV) 


1933 


2118 


5119 


5304 



We have: 

(17) m [cs] = m [cq] -m q + m s 

m[bq] = m [cq] - m c + m b 
m[ bs \ = m [bq] -m q +m s 

The numerical values for the diquark masses are given in Tab. IV. 

The last step for the mass spectrum determination is the diagonalization of the Hamil- 
tonian in Eq. (12). 

As done in [17] we label the particle states with the notation |5 q , S§\ 5 qq , J) where q 
and q represent diquark and antidiquark states respectively and S' qq is the total spin of 
the diquark-antidiquark system. In the most general case [91(72] [9394] we have: 

• two positive parity states with J p = 0+ 

(18) [0+K - |0 9l92 ,0 ?39 - 4 ;0 qq , J = 0) 

(19) |0 + )2 = |l<Jl<?2! l?394i qq , J = 0) 



• three states with J = 1 and positive parity 



(20) 


|l + >r = 


|1<I1<72> 0<J394> 


Iqq ! J = 


= 1) 


(21) 


|l + >2 = 


|0qi«2 1 1?3<34 ! 


Iqq ! J = 


= 1) 


(22) 


|l + >3 = 


|l|Jl<?2> 19394 1 


Iqq 7 J ~ 


= 1) 



We observe that if the diquark and the antidiquark have the same flavor composition, 
i.e. [<7i<Z2][<Zi<72], the charge conjugation is a symmetry of the state. Thus we define: 

• two positive parity states with J PC = ++ 

(23) 10++^ = |0 9l92 ,0 g - l9 - 2 ;0 qq ,J = 0) 

(24) |0++) 2 = |lgi92! l?i<72! qq , J = 0) 



one state with positive parity and positive charge conjugation and J = 1 

(25) |1 + + ) = (|1<?i92 ' 0«i92 i Iqq ' J = 1) + |0gi92 ) l?i«2 I Iqq ' J = 1)) 
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• two states with positive parity and negative charge conjugation 

(26) |1 + ) 1 = —^=(\l qi q 2 ,0q 1 g 2 ;l^,J=l) — \0q iq2 ,lq 1 q 2 ',l^,J=l)) 

(27) |1 + ) 2 = | lgiqr 2 ; lgi«j 2 i Iqq i ^ = 1) 

As introduced in [3] the action of the spin operators in Eq. (14) can be evaluated using 
matrix representaion for the states given above. Individual diquark spins are expressed 
by the 2x2 Pauli matrices 

(28) r° = -L„ a r = 

for spin and spin 1, respectively. The matrices T are normalized in such a way that: 

(29) Tr[(r Q ) t r /3 ] = 8 aP 
It is straightforward to write: 

|0 q , q ; qq , J = 0) = i<7 2 ® cr 2 
|lq,lq;0 q q, J = 0> = ^= (<7 2 cr l ) ® (<72ff*) 

|0q, Iq! Iqq, J = 1) = -<7 2 ® (ff2CT l ) 
|lq,0 q ; l qq , J = 1) = ^ (cr 2 (7 l ) ® <7 2 

(30) |lq, Iq; Iqq, J = 1) = ^£ l3k {<T2<T j ) ® (a 2( 7 fc ) 

The states we want to focus on are hidden charm [cq] [cq] and hidden bottom [bq] [bq] 
states, where q is a light quark u, dors. A useful way to visualize the possible tetraquark 
states is to organize them in multiplcts with defined J p ( c ) quantum numbers showing 
the mass value versus the third isospin component. All the plots shown in Fig. 1 have the 
same range in both axis in order to better visualize the mass splittings for the different 
configurations. Since the diquarks [cu] and [cd] are assumed to have the same mass, we 
will also have degenerate states. In Fig. 1 the following multiplets are shown: 

J p = + : contains 18 states. These states are taken as L = and S = states. This 
configuration requires the diquarks and the antidiquarks to have the same spin. 

Starting from the lightest masses, we find two iso-triplets, one with both diquarks 
in spin 1 and one with both diquarks with the spin 0. The corresponding two iso- 
singlets are degenerate in mass. We observe that the pink (triangle) states in the 
figure are charged, while the light blue(star) are neutral and have positive charge 
conjugation quantum number C = +1 (see Eq. (24)). 

The four red (rectangle) states shown are open strange iso-doublets: here the same 
mass value for [cu] and [cd] leads to the same mass value for the charged states 
(like [cu][cs]) and the neutral one (like [crf][cs]). 
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The two states shown in dark blue(box) are the two iso-singlets [cs] [cs] with definite 
J PC = ++ quantum numbers. 

J p = 1 + : there are 27 states. Assuming L = 0, there are three different spin configura- 
tions that give J = 1 as defined in Eq. (22). 

Starting from the lightest states, we find three iso-triplets and the correspond- 
ing iso-singlets degenerate in mass. Here again we have charged states shown in 
pink(triangle), while the neutral states are drawn in light blue(single star) if C = +1 
(see Eq. (25))and in green(double star) if C = -1 (sec Eq. (27)). 

In the middle of the plot, in red(rectangle), the three open-strange iso-doublets are 
shown, each of which has a corresponding charge state degenerate in mass. 

The three heaviest states (hidden-charm hidden-strange) are drawn in dark blue(singlc 
box) for J PC = 1 ++ , in green(double box) for J PC = 

J p = 2+ : contains 9 states. These states are taken as L = and S = 2 states. 

Starting from the lightest masses, we find one iso-triplet (both diquarks are in spin 
1) and the corresponding iso-singlet is degenerate in mass. The pink(triangle) state 
in the figure is charged, while the light blue(star) is neutral and has positive charge 
conjugation quantum number C = +1. 

The two red(rectangle) states shown are open strange iso-doublets: each one is 
double degenerate in mass since the assumption mui = mt c( {\. 

The state shown in dark blue(box) is the iso-singlct [cs] [cs] with definite J PC = 2 ++ 
quantum numbers. 

J p = + radial excitation: the same as the J p = 0+ mass spectrum is here repro- 
duced adding to each mass value a correction to obtain the second radial excita- 
tion. The only information we make use of is the mass spilitting AM = 0.657 GeV 
between the r] c (2S) and T] C (1S). Having this states the same quantum number as 
the tetraquark multiplet under investigation wc simply evaluate the second radial 
excitation adding the mass splitting AM = 0.657 GeV to each J p = + ground 
state. 

J p = 1+ radial excitation: the radial excitation splitting between the states ip(2S) 
and J/ip(lS) is used to evaluate the mass spectrum for the tetraquark multiplet 
with these quantum numbers. We evaluate numerically the mass spectrum for the 
J p = 1 + second radial excitation adding a mass splitting AM = 0.589 GeV to each 
J p = 1 + ground state. 

The same analysis can be done for the hidden-bottom tetraquark states and the 
resulting multiplcts are shown in Fig. 2. The conclusions discussed above hold for the 
[bq] [bq] system replacing charm with bottom. The only difference is that for the radial 
excitations where the only information we have to compute the radial excitation gap 
comes from the splitting AM = 0.563 GeV between the states T(25) and T(15). This 
value is used independently on J. 

The relativistic approach for L ^ 

When the two diquarks are orbitally excited spin interactions are suppressed by the 
(5-Dirac behaviour in Eq. (15). In this case we can study the mass spectra using a 
relativistic string model. 
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Figure 1.: Hidden charm multiplets with zero orbital excitation (L = 0): on left the first 
radial excitation, on right the second radial excitation. 



If orbital excitation is switched on the system can be visualized as a relativistic string 
spinning with angular velocity u connecting two masses mi, to 2 at distances n and r2 
away from the center of rotation respectively. This is the generalization of the Regge 
formula relating the masses of the hadrons M to their angular momenta L: 



(31) 



M = o-L 



Eq. (31) arises from solving the equations for a relativistic spinning string with a tension 
T = (t/2tt, terminated by the boundary condition that both ends move transversely at 
the speed of light. 

Let us consider a system composed of two masses m, connected by a string spinning 
with angular velocity u — Vijri. Each single mass experiences a relativistic central 
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Figure 2.: Hidden bottom multiplets with zero orbital excitation (L = 0): on left the first 
radial excitation, on right the second radial excitation. 



acceleration, direct along the string axis and given by: 

(32) T = mi^r^i 2 
where 7^ it the relativistic factor 

(33) 7 * ' 



The tension force T contributes to the energy of the system with 

T d 

(34) dE i =Tf i dr i 



where the relativistic ji factor arises in the laboratory reference frame. Thus the total 
energy of a spinning string is the sum of the rest energy and the rotational energy of 
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each mass m^: 

(35) E 8tring = mi7i + m 2 7 2 + - / , dv + - / di;. 

^ Jo V 1 — « w Jo vl — v 

In the same way the orbital angular momentum is evaluated. The string contribution is 

T Vi 2 

(36) dLi = Lodli = wvi dEi = —= 1 =dvi 

to V 1 — vc 

and adding the single mass terms we have the total angular momentum: 

J> run v 2 rp rur 2 y 2 

(37) L = miwr?7i + m 2 Lor\ l2 + / ~ 7 ==dv + / -==dv. 

u Jo vl-" w Jo y\ — v z 

After solving the integrals 

T 

(38) E 'string = ?™i7i + m 2 7 2 H (arcsin[wri] + arcsin[wr 2 ]), 

(39) L = miLorlfi + m 2 u;r 2 7 2 + 

— ^ ( — wri v/l - (wri) 2 + arcsin[wri] - wr 2 J\ — (wr 2 ) 2 + arcsin[u;r 2 ] N ) 
or 2 V / 

An important observation is that hadrons have not only orbital angular momentum 
but also spin angular momentum. In order to describe in more detail the hadrons mass 
spectra we can introduce spin-orbit interactions as: 

(40, EsL=k i.(^y.3 

where k is a free parameter of the model and R = n + r 2 . The total energy of the state 
will then be given by the sum of the two contributions: 

(41) E = E atring + E SL 

The free parameters of this model are: the masses mi and m 2 , the string tension 
T = a/2-K and the spin-orbit coupling k. Given this input, for each fixed J, L, S quantum 
numbers, the distances r\ and r 2 and the angular velocity w are defined. 

The free parameters are evaluated applying this relativistic string model to qq mesons 
with L ^ 0. We use as input the cc and bb mesons for the hidden charm and hidden 
bottom tetraquarks states respectively. It is surprising how the fitted values for charm 
and bottom quarks in this model are the same as in the constituent quark model, within 
a few MeV. This allows us to use in the relativistic string model the diquark masses 
already estimated in the framework of a non relativistic constituent quark model. 

In Fig. 3 the agreement between the experimental data (orange box) and the corre- 
sponding fitted mass values is shown. 

Using as input the \c iso-triplet, the h c iso-singlet and the ^i>(3770) mesons we ob- 
tain the following parameters used for the determination of the hidden- charm [cq] [icq] 
tetraquark mass spectra (L = 1): 
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Figure 3.: Standard cc and bb mesons: agreement between experimental data and theo- 
retical evaluation. 



,. nS T = 0.159 GeV 2 . . . , , 

(42) fc = 0.214 GeV- 2 ^ hldden charm 

For the hidden bottom sector we have the \b iso-triplet and the T(1D) mesons. Fit 
results give: 



T = 0.261 GeV 2 , . . . . , 
(43) fc = 0.046 GeV" 2 ^ hlddcn bottom 

It is remarkable that the fitted values for the string tension T = a/2ir in both heavy 
systems are very close to what is found by Wilczek and Selem in their phenomenological 
analysis of Regge trajectory in light mesons states [16]. 

In Fig. 4 are organized the [cq][cq] multiplets. The orbital excitation L — 1 leads to 
negative parity states. We have: 

J p = 0~ multiplet is realized only if the total spin of the tetraquark is S = 1. Since 
the model depends only on the total spin (and not on the single diquark spins) 
we have the same mass value for [cg]s = o[cg]s = iand [cq]s=i[cq]s=ista,tes. For this 
reason each shown state is a double-degenerate state. Different colors (symbols) 
refers to different flavor compositions: starting from the lightest states we have in 
light blue(star) a iso-triplet with no strange quarks, at higher mass values we find 
an iso-doublet with open-strangeness drawn in red (rectangle) and the heaviest state 
is a hidden-strange iso-singlet drawn in dark-blue(box). 

J p = 1" multiplet is realized with three different total spin angular momenta S = 
(double-degenerate states), S = 1 (double-degenerate states), S = 2 and the cor- 
responding mass splitting is shown for each different flavor composition: in light- 
blue(star) we have no-strange states, in red(rectange) we have open-strange states, 
in dark-blue(box) we have hidden-strange states. We observe that, as in the h c ~Xci 
system, the mass values increase with decreasing tetraquark total spin. 
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J p = 2~ multiplet is realized with two different total spin angular momenta: S ~ 
1 (double-degenerate states) and S = 2. The same color legend as for J p = l~ is 
used and again we observe that the mass values increase with decreasing tetraquark 
total spin. 

radial excitation multiplets with L = 1 are estimated assuming the splitting AM = 
0.360 GeV observed in the Xb multiplets for both hidden charm and hidden bottom 
tetraquark mesons. 
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Figure 4.: Hidden charm multiplets with L—l: on left the first radial excitation, on right 
the second radial excitation. 



2T.2. Production and Decays. As shortly discussed in [18] there are different decay 
modes for a given tetraquark state. For heavy states, like hidden charm and hiddem 
bottom, the most important are: 



• quark-antiquark pair production: the final state will be a two particle state, a 
tetraquark and a qq meson; 
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Figure 5.: Hidden bottom multiplets with L=l: on left the first radial excitation, on right 
the second radial excitation. 



• quark-antiquark exchange diagram: depending on which quarks exchange we have 
different two-mesons final states; 

• stretching the color string that bind the diquark to the antidiquark, a quark- 
antiquark pair is formed, letting two baryons in the decay products. 

2'2. Molecules. - Among the possible multiquark states, the main alternative to 
tetraquarks are hadron molecules, bound states of two (or more) hadrons. They were 
proposed a long time ago to describe the deuteron as a bound state of nucleons. After- 
wards, many of the resonances which do not fit in the standard qq and qqq scenario have 
been given a molecular assignment. Nevertheless, it is challenging to identify possible 
multiquark states in an environment of many broad and often overlapping conventional 
states. More reliable experimental indications of the existence of exotic hadronic can- 
didates came recently when searching for more cc mesons. The first hadron molecule 
interpretation in the charm meson sector was proposed for the tp(4OA0), observed in 
e + e~ annihilations [19, 20]. 
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The main difference between molecules and tetraquarks is that the hadrons bound in 
the molecule preserve their "atomic integrity" . This obviously leads to different predic- 
tions for spectra, production mechanisms and decay rates. Before analyzing the details 
of the dynamics we give some distinctive signatures of hadronic-pair candidates [21]: 

1. in the light sector J PC and flavor quantum numbers should be compatible with 
those of an L = hadron pair (nuclear forces which bind the two mesons together 
are short ranged so that a P-wave bound state seems unlikely). 

In the heavy mesons sector there is the possibility of P-wave bound states but this 
would imply the existence of more deeply bound S'-wave molecules [22]. 

2. binding energies should be of the order of 50 McV for the light mesons and of the 
order of 10 MeV for the heavy ones. This is because the minimum distance required 
for the hadrons to maintain separate identities is R ~ 1 fm and the binding energy is 
Eg ~ l/(2/iP 2 ) (where /i is the reduced mass). As an example KK molecules have 
H ~ 500 MeV and thus E B < 50 McV, whereas DD molecules have fj, ~ 1 GeV 
and thus Eb < 10 MeV, and finally for BB molecules fi ~ 5 GeV leading to 
E B < 4 MeV. 

3. Large branching ratios into final states containing the constituent mesons despite 
the reduced phase space available. 

4. Anomalous electromagnetic coupling with respect to the ordinary charmonium 
states. 

The spectrum of hadron bound states can be predicted once a model for the interaction 
between the hadrons has been proposed. The interaction between two hadrons varies with 
the distance: at short distance quarks inside the hadron interact with each other through 
the exchange of gluons, at long distance the exchange of mesons (tt, lu, p, ...) between 
the hadrons themselves is dominant. We now discuss separately these two dynamical 
regimes. 

Gluon exchange. Hadrons arc bound states of quarks, which interact with each other 
through electric and color charge. Thus, one can study the problem of hadron-hadron 
interactions and consequently hadron-hadron binding by looking at their constituents 
interactions. This interactions is mediated by the gluons, which are the vector boson 
particles associated with the gauge fields of the SU(3) non abelian gauge field theory 
which quantitatively describes the strong interactions. 

At first order in perturbation theory using g s (the strong coupling constant) as per- 
turbative parameter, we are able to compute the non relativistic potential associated 
with one gluon exchange diagrams. This results in a Coulombic term, to which a linear 
confining term must be added to account for the confinement in the strong coupling limit 



where is the distance between the i-th and j-th quarks or antiquarks which constitute 
the hadrons and b is a string tension parameter. The confinement is indeed phenomcno- 
logically implemented through a string model. The string is responsible for the attractive 
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force between two quarks at large distances. The energy associated to a string with two 
masses at its ends is oc br, b is the string tension. 

Relativistic corrections to the potential can be included as in the case of the hydrogen 
atom, just introducing in the calculation the SU(3) generators A a . Through a perturba- 
tive expansion in (p/m) 2 one is able to derive the fine and hyperfine structure equivalent, 
in the atomic model, to the spin-orbit and Darwin terms and contact spin-spin interac- 
tions, respectively. The two terms, in the case of strong interactions read: 



The hyperfine term can be understood as a dipole-dipole interaction: to each quark 
is associated a chromo-magnetic dipole moment which is proportional, just as in the 
electron case, to the coupling g s and to the inverse of its mass. This spin-spin term is 
weighted by a smeared delta function, usually a gaussian, which accounts for the short 
distance nature of the interaction. The interaction potential is, thus, the same as that 
used successfully to describe the conventional meson/baryon spectrum. The parameters 
of the interaction hamiltonian can be thus derived from ordinary meson spectroscopy, 
and there is generally a quite broad agreement in the literature, as summarized in Tab. V. 



Table V.: One Gluon Exchange parameters. 





m„(GeV) 


m s (GeV) 


m c (GeV) 


a c 


C (GeV) 


6(GeV) 




a (GeV) 


[23] 


0.375 


0.600 




0.748 


-0.777 


0.178 






[24] 


0.375 


0.650 




0.857 


-0.4358 


0.154 


0.840 


0.70 


[25] 


0.335 


0.550 


1.600 


0.590 




0.162 


0.590 


0.9 



Meson exchange. The existence of bound states of two mesons due to one pion ex- 
change was proposed for the first time by Tornqvist in [26]. The inspiration came from 
the deuteron, a bound state of two nucleons which interact through the long range One 
Pion Exchange potential: 



2 p-m, 

(46) V n (f) = ^ (n ■ r 2 ) [(3(*i • f) (<? 2 ■ f) - u x ■ ff 2 ) W(r) + (a, ■ <r 2 )} — - 
6 r 



with 



(47) W(r) = 1 + - JL, + — 

where 01,2 are the spins of the nucleons, and T\ t 2 their isospins. The potential contains 
a scalar term proportional to the spin-ispospin factor with a pure Yukawa potential, and 
a tensor term, which is a higher order correction. Let's discuss the deuteron case, since 
it will be useful for comparison later on in our discussion. The spin-isospin factor gives 
an indication of the potentially binding channels: 
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+9/16 (5 = 0,7 = 0) 
-3/16 (S = 0,7 = 1) 
-3/16 (S = 1,7 = 0) 
+1/16 (5 = 1,7= 1) 

The binding is expected in the (S = 0, 7 = 1) channel and in the (S = 1,7 = 0) 
channel, the deuteron one. The first possibility is ruled out when including also the 
tensor term, which instead strengthens the attraction in the second configuration adding 
a 7>-wave component. 

Tornqvist speculated on the possibility that such a potential could bind pairs of 
mesons, calling deusons these deuteron-like bound states. Since the pion is very light, 
deusons can be very large, much larger than ordinary qq mesons. For ground state mesons 
there are two possibilities: pseudoscalar-vector mesons and vector-vector mesons, since 
bound states of pscudoscalar-pseudoscalar mesons are forbidden by parity. When consid- 
ering 5- wave bound states, the relative strength and sign of the potential for different spin 
and isospin channels is given by the so-called relative binding number, firstly introduced 
by Tornqvist in [26], which is the analog of (48) for meson pairs. Thus the relative bind- 
ing number gives an indication of the binding channels. If one allows for higher angular 
momentum the spectroscopy becomes more complicated and more detailed calculation 
are needed. Moreover in [22, 27] the possibility for S-wave pions is considered. 

Manohar and Wise studied in [28] the interaction between two heavy-light mesons 
Qq and considered the possibility that a BB bound state could exist. Their derivation 
of the potential proceeds from the Heavy Quark Effective Theory formalism combined 
with the chiral Lagrangian approach to describe the light pseudoscalar mesons. Their 
potential agrees with (46) . 

The same calculation scheme is applied by Tornqvist in [29] to meson-antimeson 
bound states (remind that the meson-antimeson interaction potential is the opposite of 
the meson- meson one). The only parameter of the model, besides the mesons masses, 
is the pion-meson-meson coupling g. Its value can be deduced from the ttN coupling 
exploiting the relation 



(48) (n • f 2 ) (ai • <r 2 ) 



S(S + 1) 



(49) ^ = ^%™ 



9%N_ = ^S^_2 

4tt 9 f 2 



where / is the decay constant of the pion. This relation arises when one uses chiral 
perturbation theory to describe the pseudoscalar octet mesons. In this way one can 
derive the effective quark-pion interaction 

(50) C lnt = jq^fqd^TT 

Treating the quark fields as constituent quarks and using the non relativist quark model 
to derive hadronic matrix elements gives the relation (49) between the irN coupling and 
the ttDD* one. The measured value g\ N ~ 1 leads to g 2 ~ 0.6. A striking confirmation 
of this estimate is the prediction of the D* width of ^70 KeV which was later confirmed. 
The potential is computed for each of these states paying attention to the fact that there 
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will be some coupled channels (states with same J and S but different L will mix with 
each other). 

However, because of the singular nature of the tensor part of the potential, i. e. the 
first term in (46), a regularization procedure is needed. The most natural method is to 
introduce a form factor at the 7r7Vvertex, which gives to the pion source a spherical ex- 
tension with radius R ~ 1/A, with A an ultraviolet cutoff. Even if the phenomenological 
knowledge of the cutoff A is rather poor, it may be fixed by comparison with nuclear 
physics. From NN interactions A must be in the range 0.8-1.5 GeV, while to reproduce 
the deuteron binding energy one needs A ~0.8 GeV. The value employed in [29] and in 
the literature in general is A ~1.2 GeV, which seems appropriate for D mesons. It is 
crucial here to emphasize that the existence or otherwise of meson-meson bound states 
depends strongly on the value chosen for A [26, 24]. 

2'2.1. Spectra. We briefly review the result for the spectra of possible molecular states 
as obtained in the two approaches. 

Gluon exchange. This interaction scheme has been used to study a number of systems 
with different approaches from case to case. 

Weinstein and Isgur, by using the potential terms of Eq. (44,45), showed that /o(975) 
and ao(980) could be interpreted as a KK molecule, [30, 31, 23]. They find the full 
four-quark wavefunction of the bound states using a variational Gaussian basis. Then 
they invert the Schroedinger equation to obtain an effective potential and integrate this 
potential to obtain the observed phase shifts in nir scattering for ao and / . This model 
predicts that in general the qqqq (with q — u,d, s) ground states are two unbound mesons 
except for a KK molecule. 

Another case in which the one gluon exchange has been successfully exploited is that 
of A (1420) as a K* K bound state [32, 33]. 

Furthermore in a series of papers by Swanson and Barnes [24, 34, 35] meson-meson 
interactions have been obtained extracting an effective potential between the two mesons 
from the Born order scattering amplitude computed from Eq. (44,45). Such interactions 
necessarily involve the exchange of quarks between mesons. This is because the appli- 
cation of the one gluon exchange hamiltonian induces transitions among different color 
configurations. For example (lj®lj|A"A°*|8j®8j) ^ 0. When searching for molecules one 
wants to deal always with two color-singlet objects and thus it is necessary to exchange 
quarks or antiquarks between the two mesons to obtain again a 1 ® 1 configuration. 
Nevertheless this consideration has a conceptual drawback, as explained in [24], since if 
quark exchange must occur the range of the effective potentials is limited by the mesonic 
radii to roughly 1 fm. Thus a two-meson bound state obtained from one gluon exchange 
has a spatial extension which cannot exceed 1 fm: at this point the difference between 
hadronic molecule and tetraquark state seems to be just a matter of language, the only 
difference between the two being the way in which color is saturated. 

Meson exchange. Bound states of mesons due to one pion exchange have been studied 
in [28, 29] and many others. 

In [28] the authors compute the binding energy of an eigenstate of the potential 
using a variational calculation. They assume that the short range potential, which is 
dominated by vector meson exchange, cannot be repulsive and find that the BB* in the 
(I = 0, S = 0) channel leads to a binding energy of 8.3 MeV (for g 2 — 0.5). The authors 
state that the value of the binding energy is sensitive to the precise value of g 2 and 
to the value chosen for the potential at short distances. They neglect the contribution 
of heavier mesons exchange, arguing that they are less important than in the nucleon- 
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nucleon case. On the other hand no DD or DB bound state is found because of the 
positive contribution from the kinetic energy which overwhelms the attraction due to the 
potential. 

In [29] the method used to find bound states is to solve numerically the Schroedingcr 
equation with the one pion exchange potential of Eq. (46). The author considers flavor 
non exotic states in the heavy- heavy (Tab. VI) and light-light sectors (Tab. VII), and the 
flavor exotic heavy- light and light-light states (Tab. VIII). The general pattern is that 
for pseudoscalar- vector mesons the binding is more likely in J PC = h , I ++ , while for 
vector- vector ones in J PC = h , ++ , 1 H , 2 ++ . In both cases the isospin configuration 
is I = 0. 



Table VI.: Flavor non exotic bound states in the heavy sector found in [29]. The 
£)(*)£)(*) are predicted to be almost at threshold, compatible with zero binding energy, 
while for the £?(*).£?(*) molecule the typical binding energy is of the order of 50 MeV. 



Constituents 


M (MeV) 


J PC 


A(GeV) 


DD* 


~ 3870 


0"+ 


> 1.5 




~ 3870 


1++ 


1.2 


D*D* 


~ 4015 


0"+ 


1.5 




~ 4015 


0++ 


1.2 




~ 4015 


1+- 


1.3 




~ 4015 


2++ 


1.2 


BB* 


~ 10545 


o-+ 


1.2 




~ 10562 


1++ 


1.2 


B*B* 


~ 10590 


0"+ 


1.2 




~ 10582 


0++ 


1.2 




~ 10608 


1+- 


1.2 




~ 10602 


2++ 


1.2 



In [36] Swanson obtained the same results for the vector-vector heavy-heavy mesons 
bound states, except for the fact that he finds only one D*D* bound state with J PC = 
0++. 

In Tab. IX the possible molecular assignments for exotic meson candidates are sum- 
marized. 

2 '2. 2. Production and decays. The production of molecular bound states proceeds 
necessarily through the simultaneous production of its constituent mesons in a suitable 
relative momentum configuration. On general grounds it is reasonable to assume that 
the relative momentum between the two mesons which are candidates for the binding 
is related to the binding energy of the molecule itself. This relation can be deduced 
using the uncertainty (or minimal uncertainty) principle: (Ax)(Ap) ~ h. The binding 
energy gives an estimate of the size of the bound state Eb = h 2 /2fi(Ax) 2 and thus the 
average momentum spread is (Ap) ~ ^/2/iEb- Each of the molecules predicted nearly at 
threshold, thus with a tiny binding energy, will allow for very small momentum spreads. 
A rather special case is the A" (3872) which will be treated in detail later. 

The decay modes of a meson-meson bound state can be divided into two classes: 
long-distance decay modes and short- distance decay modes. 
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Table VII.: Flavor non exotic bound states in the light sector found in [29]. "waa" 
stands for "weak attraction added" , which means that the binding is obtained if a weak 
short-range interaction potential (~ 7j e~( r / r °' ) is added to (46). 



Constituents 


Threshold (MeV) 


Candidates 




A(GeV) 


BE expected (MeV) 


KK* 


1394 


77(1440) 


o-+ 


waa 


~ 50 






/i(1420) 


1++ 


waa 


~ 30 


K*K* 


1792 


77(1760) 


o-+ 


waa 


~ 30 






/o(1710) 


0++ 


waa 


~ 80 








1+- 


waa 








/ 2 (1720) 


2++ 


? 


~ 70 


(pp + 


1550 - 1566 


77(1490) 


o-+ 




~ 70 


(pp - uu)/V2 


1550 - 1566 


77(1515) 


0++ 




~ 40 


(pp + uju)/\/2 


1550 - 1566 


/ 2 (1520) 


2++ 




~ 40 


{K*p- K*u)js/2 


1671 - 1679 




0++ 







The former class consists of the decaymodes of the constituent mesons. The partial 
decay widths of the molecular state in these modes are related to those of its constituent 
mesons. If the binding energy is very small, nearly zero, the partial decay widths of the 
molecule in these modes will be almost equal to the ones of its constituents, whereas for 
deeply bound states one expects large deviation from the free meson widths. In particular 
a large binding energy tends to stabilize the meson when it is bound. 

The latter class of decay modes, that of short-distance decay, is a manifestation of 
the existence of some inelastic channels in the meson-meson scattering mechanism. The 
associated partial decay widths are proportional to the probability that the two mesons 
come together at a point, which is given by the square modulus of the bound state wave 
function at the origin |i/>(0)| 2 . 

2'3. Hybrids. - Both the molecular and tetraquark hadrons are built with quarks 
and antiquarks. However the QCD Lagrangian contains also the gluons, which can 
act as dynamical degrees of freedom besides being the particles which mediate strong 
interactions. One can indeed suppose the existence of gluonic hadrons, bound states 
of gluons and quarks. There are two kinds of gluonic hadrons: the glueballs, which are 
bound states of only gluons, and the hybrids which are qqg bound state, that is a qq state 



Table VIII.: Flavor exotic bound states found in [29]. They are expected in the isoscalar 
CP and 1 + channels and in the isovector 1~ channel. 

Constituents Threshold (MeV) 

KK* 1394 

DD* 3871 

BB* 10605 

KD*+K*D 2505 - 2761 

KB*+K*B 5823- 6176 

DB*+D*B 7190- 7286 
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with a gluonic excitation. This is is not surprising from the point of view of color, since 
qq € 8 © 1 and g £ 8 and one can pick a singlet component from the 8 <g> 8 configuration. 
In this section we will focus on hybrid mesons. 

The existence of hybrid mesons was suggested in 1976 by Jaffe and Johnson [85] and 
Vainsthein and Okun [86]. Hybrids have been studied using different approaches: (i) the 
MIT bag model, (ii) an adiabatic heavy-quark bag model, (m) constituent gluon models 
and (iv) heavy quark lattice gauge theory, (v) the flux tube model, The average mass 
obtained for the lightest hybrid with light quarks is about 1.5-2 GeV. Hybrids can exhibit 
exotic J PC quantum numbers, and thus could be easily identified experimentally. 

The MIT [87, 88] bag model predicts the existence of a lightest hybrid mesons mulit- 
plct at ~1.5 GeV and the presence of an exotic 1 h state in this multiplct. The exotic 
J PC quantum numbers are due to the boundary conditions in the bag. 

For the heavy quarks a spherical bag would be quite unrealistic, and thus an adiabatic 
bag model was introduced by Hasenfratz, Horgan, Kuti and Richard in [89]. In this 
model the bag was allowed to deform in the presence of a fixed QQ source. The resulting 
potential was used in a Schroedinger equation to compute the mass of the hybrids. The 
lightest hybrids for cc was found at ~3.9 GeV and at ~10.5 GeV for bb. For a recent 
result on adiabatic potentials in QCD string models see [90]. 

Constituent gluon models treat the gluon as the constituent quark model treat the 
quarks. These models have been introduce by Horn and Mandula in [91] and later 
develped by Tanimoto, Iddir et al. [92, 93, 94] and Ishida et al. [95, 96]. The gluon 
has a fixed orbital angular momentum relatively to the qq pair, usually called l g , and 
the cc is in a defined orbital configuration l qq and spin configuration s qq . The quantum 
number of such a bound state are somewhat different from the ones predicted in other 
model: P = ( — l)'a+^<? and C = (—l) l n+ s n+ 1 . The lightest hybrid states in this model 
has l g = and thus non exotic quantum numbers such as 1 are obtained using p-wave 
qq states with s qq = 1, while exotic 1 h state has s qq = 0. 

Lattice QCD is supposed to give the most reliable predictions for absolute hybrid 
masses. In the heavy quark lattice QCD in which the QQ pair is kept fixed while the 
gluonic degrees of freedom are allowed to be excited the lightest charmonium hybrids 
was predicted in [97] to have a mass of 4.2 GeV. for cc and 10.81 GeV for bb. 

Finally, the flux-tube model is the most widely used model for hybrids. In lattice QCD 
two separated color sources are confined by approximately cylindrical regions of color 
fields if they are sufficiently far apart. The flux tube model describes this confinement in 
a simple dynamical way, approximating the confining region between quarks by a string 
of massive points. This approach is motivated by the strong coupling expansion of lattice 
QCD. Since in a lattice gauge theory the flux line can be expanded only in transverse 
directions, in the flux tube models one allows only for transverse spatial fluctuations of 
the massive point positions. In the first studies with this model an adiabatic separation 
of the quark and gluon degrees of freedom was carried on, exploiting the fast dynamical 
response of the flux tube degrees of freedom with respect to the heavy quarks time 
scales. This separation allows to fix the qq separation at some value R and compute 
the eigenenergy of the system in some fixed configuration of the flux tube: E\(R). The 
ground state A = gives the ordinary meson spectrum. Hybrids are obtained for A > 
and can be studied using the excited potential E\(R). The lightest hybrid state is the one 
in which the string has a single orbital excitation about the qq axis. In initial models the 
adiabatic potentials were determined in the approximation of small fluctuations relatively 
to the qq axis. This approximation was removed by Barnes, Close and Swanson in [98]. 

In the charmonium family hybrids are predicted in the mass region 4.3 GeV, with an 
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estimated uncertainty of 100/200 MeV. As for the bottom sector hybrids are predicted 
in the region 10.7-11.0 GeV. 

While the masses of hybrid mesons are computable in all the models listed above, and 
in particular in Lattice QCD, the decay dynamics is more difficult to study. 

The only model which offers a description of the decay dynamics is the flux-tube 
model. Indeed in this context the decay occur when the flux-tube breaks at any point 
along its length, producing in the process a qq pair in a relative J PC = ++ state. A 
similar model has been applied to the ordinary QQ mesons, since, as we stated before, 
the flux-tube model in its ground state describe ordinary mesons. The distance from the 
QQ axis at which the light pair is created is controlled by the transverse distribution 
of the flux-tube. This distribution varies when going from the non-excited flux-tube to 
the first excited flux-tube configuration. Exploiting the empirical success of this model 
in describing the ordinary mesons decay dynamics Close and Page in [99] derived the 
decay pattern for hybrids. They found that in a two meson decays the unit of orbital 
angular momentum of the incoming hybrid around the QQ axis is exactly absorbed by 
the component of the angular momentum of one of the two outgoing mesons along this 
axis. In [99] they treated explicitly the light flavor case, but a generalization to hybrid 
charmonia is straightforward. The final state should be in this case £)(*,**)£)*,** ^ where 
D** indicates D-meson which are formed from p-wave cq (q = u,d) pairs. However, 
since the masses predicted in the flux-tube model are about ~4.3 GeV, i.e. below the 
DD** threshold, it is possible that this decay is kincmatically forbidden giving a rather 
narrow resonance decaying in charmonium and light hadrons. These modes offer a clear 
experimental signature and furthermore should have large branching fractions if the total 
width is sufficiently small. 

2'4. Hadrocharmonium. - Recently [100] a new interpretation has been proposed for 
the states with J PC = 1~ in the region 4.2 - 4.6 GeV, namely Y (4260)-F(4350)- 
y(4660), and for the only charged state observed at 4.43 GeV, namely the Z + (4430). 

These states show some common characteristics: they decay prominently in only one 
of the two charmonium states J/ip and tp(2S) and furthermore the decay into open- 
charm mesons is highly suppressed. These common features have been interpreted as the 
indication of an hadronic structure in which a standard charmonium state is stuck into a 
light hadrons. This picture is inspired by the much discussed case of charmonium states 
bound inside a nucleus. This charmonium state embedded inside light hadronic matter 
is referred to as hadro-charmonium or hadro-quarkonium in general. The light hadronic 
matter act as a spatial extended environment in which the more compact J/i[> or ip{2S) 
moves. This picture is at least able to explain why the decay into J/ip or ip(2S) is favored 
or suppressed, depending on which charmonium state is stuck inside the hadron. 

The reason why the cc state interacts with the light-hadronic stuff although being 
neutral with respect to color charge, is that it possesses a chromo-electric polarizability. 
Thus its chromo-electric dipole moment interacts with the chromo-electric field generated 
by the light hadronic matter. This interaction can be treated with the multipole expan- 
sion in QCD used for the charmonium binding inside nuclei. The chromo-electric dipole 
moment is proportional to the chromo-electric field E a through the chromo-polarizability 
a, resulting in an effective interaction hamiltonian of this form: 



(51) 



H e ff — --«E" • E a 
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The chromo-polarizability a can be deduced from the decay ip — > J/ifm + Ti~ . The 
hadronic transition in quarkonium systems arise from the interaction of a quark or an 
antiquark with the gluons, which can the materialized in light hadrons, such as 7r or r\. 
This kind of interactions can be treated within the multipole expansion, as it happens 
for the interaction of heavy quarks and antiquarks with photons [101]. The radiative 
transition terms are proportional to the electric dipole moment and to the magnetic 
dipole moment associated to the heavy quarks, the so-called El and Mi transitions: 



H E1 = -e c e(r-E) 

h m1 = -^ c (a-e 



where A = o\ — a^, B\ and <?2 being the spin of the quark and the antiquark respsectively. 
In the same way one can treat the interaction of a chromo-electric dipole moment and a 
chromo-magnetic moment with the chromo-electric and chromo-magnetic fields: 



H E 1 = -\C (r-E a ) 
(53) 2 

h u1 = -^? {?■*>) 

where £ a = t\ — 1% is the difference between the color generators acting on the quark and 
antiquark, and E a and B a are the chromo-electric and chormo-magnetic component of 
the gluon strength tensor. The two pion transition ip' — > J/ipir + Tr~ is generated by two 
insertions of the operator Hei in Eq. (53): 



1 



(54) H eff = -^ 2) E?E? 

where 



(55) a™ = -(lS^nQr^S) 

Q being the two point Green function of the heavy quark pair in a color octet configura- 

( 12) 

tion. In the leading non relativistic order for transitions in 5-wave a\- actually reduces 
to a scalar a' 12 -* , wich is measured to be a' 12 -* ~ 2 GeV -3 . On the other hand the average 
value of the product of chromo-electric fields over the light hadron X can be estimated 
using the conformal anomaly relation in QCD: 



(56) (X\-E a -E a \X)>—M x 

In this way one is able to estimate the strength of the interaction between the light 
hadronic matter and the quarkonium system bound inside it. The possibility that such 
a bound state exists depends on the relation between the mass Mx and the spatial 
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extension of the light hadron. In particular in [102] it has been shown that a quarkonium 
state does form a bound state inside a sufficient highly excited light hadron. Furthermore 
the authors in [102] find that for this kind of bound state the decay into open heavy 
falvour mesons is suppressed exponentially as exp (— ^/Xqcd / Mq) , which would explain 
the non observation of the decay of the Y resonances into final states with pairs of 
charmed mesons. 

2'5. Pentaquarks and hexaquarks. - Beyond the tetraquark structure some other mul- 
tiquarks hadrons have been considered in the literature, such as baryon states built of 
4<7 and a q, namely pentaquarks, or mesons constituted by 3q and 3q, the so-called 
hexaquarks. 

The search for pentaquarks was originally triggered by the observation of a narrow 
exotic resonance + [103] and has motivated the study of uudds state [104]. The sub- 
sequent disproof of the observation of the 9 + in photoproduction has not been taken as 
the last word. Studies on pentaquarks are ongoing. 

In particular the spectrum of the pentaquark [105] and hexaquark [106] states has been 
computed assuming that the main role is played by the chromo-magnetic interactions. 
In [105] both positive and negative parity states have been considered including also the 
splitting due to the SU(3)f breaking. In [106] the same authors have extracted the 
chromo-magnetic interaction parameters from the supposed tetraquark structures, such 
as do, /o, cr and A (3872) and use this parameters to compute the masses of qqqqqq states. 
Hexaquarks should show a strong affinity to the baryon-antibaryon decay mode. 

3. Experimental Primer 

Hadronic spectroscopy has experienced a renaissance in the last few years, thanks to 
the opportunities provided by several experiments that allowed to improve our knowledge 
of standard hadrons and discover a rich zoology of exotica. In this section we briefly 
introduce the main apparatuses and methods that arc behind the experimental results 
described in the paper. 

3T. Experiments. - In the last two decades, three e+e~ colliders, CESR [107] at 
LEPP (Cornell, USA), PEP-II [108] at SLAC (Stanford, USA) and KEK-B [109] at 
KEK (Tsukuba, Japan), have been mainly operated at a center of mass energy of 10.58 
GcV, corresponding to the mass of the T(45) bb resonance, with the main purpose of 
discovering CP violation in the B meson sector and make precision studies of CKM 
physics. 

Interesting flavour physics studies have also been performed at similar e + e~ facilities 
running around the threshold for the production of tau and charm, between 3.7 and 5.0 
GeV: BEPC at IHEP (Beijing, China) and CESR-c [110] at LEPP. 

Meanwhile the Tevatron pp collider at Fcrmilab (USA) investigated the high energy 
frontier, running at a center of mass energy of 1.96 TeV. 

All these facilities hosted experiments that demonstrated to be a perfect place for 
studying standard and exotic hadronic spectroscopy, through different production and 
decay mechanisms. 

3T.1. B-Factories. CESR, PEP-II and KEK-B, commonly named B-Factories for the 
large production cross section of B meson pairs at the T(45) resonance, operated in the 
last decades with a peak luminosity of 1.2 x 10 33 , 12 x 10 33 and 21 x 10 33 cm~ 2 s _1 , 
respectively. CESR, after collecting 15.5 fb -1 of integrated luminosity, ceased operations 
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as a £?-Factory in 1999, when PEP-II and KEK-B started providing collisions. PEP-II 
was operated until April 2008 and delivered 553 fb -1 , while KEK-B is still running and 
exceeded lab -1 . In the case of PEP-II and KEK-B, the electron and positron beams 
collide with asymmetric energies and the T(4S) resonance is produced with a Lorentz 
boost of /?7 = 0.53 and 0.46 respectively. PEP-II and KEK-B have been also run at 
the energies of the T(2S), T(3S') and T(55) resonances, in order to study the decays of 
these states. The relevant parameters of the three colliders are quoted in Tab. X. 

Three general purpose detectors were installed at these facilities: the CLEO detec- 
tor [111, 112, 113] at CESR (in three different configurations, CLEO-II/II.V/III), the 
BaBar detector [114] at PEP-II and the Belle detector [115] at KEK-B. 

The design of the the three detectors is quite similar. All of them are provided with 
a multi-layer vertex tracker of double-sided silicon strip detectors and a drift chamber, 
operating in a 1.5 T magnetic field . The silicon detectors allow a vertex resolution at 
the level of 100 ^m and are also used for standalone tracking of low momentum particles 
(Pt < 100 McV). The drift chambers are used for the tracking of higher momentum tracks 
and also provide dE/dx measurements for particle identification (PID) below 700 MeV 
and above 1.5 GeV. Momentum resolutions below the percent level are reached by these 
tracking systems. 

Different technologies have been used for PID of charged hadrons above 700 McV. 
A ring imaging Cerenkov counter (RICH) is used by CLEO. LiF is used as the radiator 
and photons are collected by a CH4/TEA detector. A aerogel-based RICH was adopted 
by Belle. Since the Cerenkov threshold for pions in the aerogel is 1.5 GeV, a time-of- 
flight (TOF) detector with 95 ps resolution is also used for K/ir separation below this 
threshold. Finally, BaBar developed a detector of internally reflected Cerenkov light 
(DIRC), where the light produced in quartz bars propagates by internal reflection along 
the bar itself and is finally detected by 11000 PMTs surrounding a water box at the 
end-cap of the detector. 

The detection and energy measurement for photons and electrons is provided in all 
three experiments by a CsI(Tl) electromagnetic calorimeter, with excellent energy reso- 
lutions, ranging from 2 to 6% in a wide energy range (from 20 MeV to 9 GeV). 

The flux return of the magnetic field was finally instrumented for the detection and 
identification of muons, by using plastic streamer counters in CLEO and resistive plate 
chambers (RPC) in BaBar and Belle. 

3T.2. Tevatron. The Tevatron pp collider at FNAL started its Run II in 2002 provid- 
ing collisions at y/s = 1.96 TeV with a typical luminosity of 10 32 cm _2 s _1 and a peak 
luminosity exceeding 3 x 10 32 cm~ 2 s^ 1 . A total integrated luminosity of 6.5 fb~ have 
been delivered up to now. 

Two general purpose detectors, CDF-II [116] and DO [117], were operated at the 
Tevatron during the Run II. Both detectors are provided with a silicon vertex tracker. 
Seven and four layers of double sided silicon strips are used by CDF and DO, respectively. 
An additional innermost layer is also included in both experiments to improve the impact 
parameter resolution, down to 40 /im in CDF. The tracking system is completed in CDF 
by a drift chamber for the reconstruction of tracks with p t > 400 MeV. In DO, a tracker 
composed by scintillating fibers read by visible light photon counters. The tracking 
systems are operated inside a solenoidal field of 1.4 T in CDF and 2 T in DO. The CDF 
tracker reached a transverse momentum resolution of <j Pt /pt ~ 0.001 ■ px [GeV] and 
allows a K/tt separation on 1.4(7 for px > 2.0 GeV thanks to the dE/dx measurement 
in the drift chamber. PID for kaons and pions below 1.5 GeV is guaranteed in CDF by 
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a TOF device with 110 ps resolution. 

Energy measurement for electrons, photons and hadrons is provided by calorimeters. 
In CDF, a segmented electromagnetic calorimeter is used, with alternate sheets of plastic 
scintillators and lead. In the hadronic calorimeter, iron is used instead of lead. Reso- 
lutions of 13.5% /y/E © 2% and 50%/^ © 3% are obtained in the central part of the 
electromagnetic and hadronic calorimeter, respectively. In DO, a LAr sampling calorime- 
ter is used, divided in three sections: a finely segmented electromagnetic section, a fine 
hadronic section and a coarse hadronic section. Uranium is used as absorber in the first 
two sections, copper and stainless steel in the third one. Energy measurement for elec- 
trons is complemented by a central preshower detector instrumented with scintillating 
fibers. An energy resolution of o E jE = 15%/VE + 0.3% and o E jE = 45%/sqrtE + 4% 
is obtained for electrons and pions respectively. 

Both CDF and DO are finally provided with muon detectors composed by tracking 
systems and scintillator trigger counters. In CDF, wire chambers are used for muon 
tracking. In DO, drift tubes are adopted and two out of three detector layers are located 
inside a 1.8 T toroidal magnetic field. 

3'1.3. t— charm Factories. Lepton colliders running slightly above 3.5 GeV allow to 
produce large quantities of tau and charm hadrons, providing an ideal environment for 
studying t physics and charm spectroscopy. The BEPC collider at IHEP covers the 
energy range between 3.7 and 5.0 GeV, with a peak luminosity of 12.6 x 10 30 cm _2 s _1 , 
with recent upgrades aimed to reach 10 33 cm _2 s~ 1 and host the BES-II experiment [118, 
119]. CESR-c is the r-charm version of CESR and covers the energy range between 3.97 
and 4.26 GeV with a peak luminosity of 76 x 10 30 cm _2 s _1 . A modified version of CLEO, 
called CLEO-c [110], is operated at this facility. Other parameters of the two colliders 
are provided in Tab. XI. 

At BEPC, the BES-II detector was composed by a straw tube vertex detector and a 
drift chamber for tracking, a scintillating TOF device with 180 ps resolution for PID and 
a sampling electromagnetic calorimeter with streamer tubes and lead absorber. These 
detector were installed inside a 0.4 T solenoidal magnetic field, whose yoke was instru- 
mented with proportional tubes for muon detection. The upgraded version of the detec- 
tor, BES-III, has a drift chamber for tracking, a better TOF with 100 ps resolution, a 
CsI(Tl) calorimeter and RPCs for muon identification. The magnetic field is increased to 
1 T. BES-II collected 58 million J/ip events and 14 million ip(2S) events, while BES-III 
started operations in 2008 and so far has collected 10pb _1 at the ip(2S). 

The CESR-c detector, CLEO-c, derives from CLEO-II, adapted to the lower energy. 
The magnetic field was reduced to 1 T and and the silicon detector was replaced by 
an inner drift chamber. CLEO-c collected 27 million of ip(2S) events, 818pb _1 at the 
V>(3770) and 602 pb -1 at a center of mass energy of 4.17 GeV, in order to study D S D S 
events. 

3'2. Analysis methods. - In this section, we briefly discuss the fundamental analysis 
approaches that are adopted in the different experiments to search for new resonances 
and investigate their properties, like masses, widths and quantum numbers. 

3'2.1. Search and measurement of mass and width. The experimental method adopted 
to look for hadron resonances depend on the production mechanism considered for the 
exotic particle (Y). The choice is among direct production, e + e~ — > Y, the initial state 
radiation (ISR) production e + e~ — > Y^jsr, 77 production e + e~ — > 77e + e~ — > Ye + e~ , 
and B decays, in particular B — s- YK^*K 
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Direct production can only be studied by scanning the center-of-mass energy of the 
accelerator (\/(s)) and requires a dedicated plan of the machine. After measuring the 
cross-section (a) cither of the inclusive hadronic production or of exclusive final states 
as a function of the center-of-mass energy, a new state would appear as a resonance in 
(j{\f\s)). Scans dedicated to the exotic charmonia were performed by CLEO-c, while 
both Belle and BaBar investigated the regions above the T(4S I ) for exotic bottomonium 
exploiting energy scans. 

In case of ISR production at e+e" colliders, it is possible, in principle, to look for 
an almost monochromatic peak in the energy distribution of the ISR photon. This fully 
inclusive method has been used in the successful search for the bottomonium ground 
state rib at BaBar [120], but it can be applied only in very special cases. In fact, there is 
a typical photon detection threshold of 30 MeV that prevents the use of this method for 
a mass near to yfs, there is a very large background for low photon energies, and photons 
are preferentially emitted along the beam axis, so that most of them escape the detection. 
Hence, most of the searches are performed by looking for fully reconstructed final states, 
requiring no additional charged particles in the event and without requiring the detection 
of the additional photon. In both inclusive and exclusive searches, the mass and width of 
the resonance can be measured. In the first case, they can be inferred from the position 
and width of the photon energy peak, and the resolution on the width is limited by 
the photon energy resolution. In the second case, the invariant mass measurement is 
available from the reconstruction of the 4-momenta of the final state particles, and the 
resonance line shape is directly observed. 

In the study of the line shape, the decay amplitude is usually described by a Brcit- 
Wigner function: 

\ m r tot (m ) 
(57) A[m) oc 



to 2 - m\ + ir tot (m)m ' 



that is possibly summed to a non resonant contribution to get the total amplitude and 
finally extract the cross section: 

(58) a(m) oc <5> } PS (m)\A{m)\ 2 , 

where $p S (m) is a phase space factor depending on the final state (for 2-body decays, 
$>p S (m) — \p(m)/p(mo)] af , where p(m) is the momentum of the decay products in their 
center of mass frame and a / depends on the final state) . Background models depend on 
the specific decay channel and will be described, where needed, in the following sections. 

Different approaches are used in the search for resonances produced in the decay 
of other particles. At the .B-Factories, standard and exotic hadrons can be produced 
in the B — > KX decays. Also in this case an inclusive or exclusive approach can be 
adopted. In the inclusive approach, one of the two -B's produced in the T(45) decay is 
fully reconstructed in a hadronic final state, and it allows to infer the 4-momcntum of the 
other B. This information, along with the 4- momentum of the kaon, allows to reconstruct 
the X mass without considering its decay products. In the exclusive approach, resonances 
are searched for in specific decay channels and the KX pair is required to be consistent 
with the hypothesis of a B meson produced from the T(45) decay. Two variables are 
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used: 



(59) 
(60) 



m ES = ^(ElY-\v* B 
AE = E* B - El , 



where is the beam energy and (E B ,p B ) is the B 4-momentum, calculated in the 
T(4S) rest frame. The "beam-energy substituted" mass rriES IS also indicated with m bc 
in Belle's papers. From the 4-momenta of the X decay products one can reconstruct the 
X invariant mass and study the line shape. 

At hadronic colliders like Tevatron, new resonances are searched for in exclusive final 
states. In particular, final states containing a J/ip decaying to [i + fx~ are searched among 
events collected with a di-muon trigger. A very large statistics is available in this case, 
allowing a precise measurement of the mass (the precision on the width is limited by the 
detector resolution) and the realization of angular analyses for the determination of the 
J PC quantum numbers. 

At the r-charm factories, the best performances can be obtained for the study of 1~~ 
states, produced in large amount by setting the center of mass energy at the mass of the 
resonance, so that the most rare decay channels can be also studied. 

3'2.2. Measurement of J PC . In some cases, the determination of the J PC quantum 
numbers can follow from the production process or the decay channels. For instance, 
resonances produced in e + e~— annihilation (with or without ISR) through a virtual 
photon can only have 1 quantum numbers. Conversely, states produced in two-photon 
events e + e~ — > e + e~-f*l* ~~ e+e~A can only have C = + and there is a set of selection 
rules for total angular momentum and parity, known as "Yang's theorem" [121], which 
hold also for resonances decaying to 77. 

In other cases, an angular analysis and the study of the kinematic spectra is needed 
in order to extract information about the quantum numbers. Considering for example 
the J/ip(fi + n~)n + n~ decay, the ir + n~ pair can be treated as a single body, so that the 
angular dependence is the one of a two body decay. For fixed helicities, it is described by 
theWi gner function D"^' ^. 2 : where Ji and are the spin and helicity of the decaying 
particle, and A^i and A ij2 are the helicities of the decay products. This dependence is 
reflected by the angular distributions of the final state particles (two muons and two 
pions in this example), usually described in terms of a few "helicity angles" (given a 
decay chain X — > A + Y with A — > B + C, the helicity angle 6a is defined as the angle 
between the B or C direction in the A rest frame and the A direction in the X rest frame) . 
Additional information can also be extracted from the invariant mass distribution of the 
7r + 7r~ pair. These analysis require large data sets and are more efficiently performed at 
the hadron colliders than the B-Factories. An example can be found in [122]. 

4. Charmonium 

The presence of heavy quarks allows very good predictions for the masses of the regular 
Charmonium states. The heavy quark inside these bound states has in fact low enough 
energy that the corresponding spectroscopy is close to the non-relativistic interpretations 
of atoms. The quantum numbers that are more appropriate to characterize a state 
are therefore, in decreasing order of energy splitting among different cigenstates, the 
radial excitation (n), the spatial angular momentum L, the spin S and the total angular 
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momentum J. Given this set of quantum numbers, the parity and charge conjugation of 
the states are derived by P = (-1)( L+1 ) and C = (-1)( L + 5 ). Fig. 6 shows the mass and 
quantum number assignments of the well established charmonium states. 




Figure 6.: Charmonium states with L < 2. The theory predictions are according to the 
potential models described in Ref. [2j. 



All the predicted states below open charm threshold have been observed with very 
good agreement with theory predictions. The Charmonium system is therefore an ideal 
environment to search for exotic states which deviate from regular spectroscopy. 

We will first review the experimental observations, reporting both the final states 
where the states have been observed and those where they have not. Since the easiest 
quantum number to attribute is C, being determined uniquely either by the production 
or the final state (see Sec. 3'2.2), we will first examine the C = + states (Sec. 41, 4'2, 
and 4'3) and then end with the J PC = 1 states (Sec. 4'4). Next we will discuss the 
evidences for charged states that play a crucial role in this field. 

In the second part of this section we will report the possible interpretations of these 
states (Sec. 4'6) and conclude by summarizing the observations and the open issues(Sec. 47). 

41. The X(3872). - The X (3872) was the first state that was found not to fit charmo- 
nium spectroscopy. It was initially observed by the Belle experiment in B — > Xifdecays 
and decaying into J/ipir + ir~ [123] and subsequently confirmed both in B decays [124] 
and in inclusive pp production [125, 126] (see the Fig. 7). This is by far the state for 
which the most information is available, and it will be here reviewed by topic: quantum 
numbers, mass and width, and production and decay. 

4' 1.1. Quantum numbers. The exotic nature of this state was initially signalled by the 
narrowness of its width (r x(3872 ) = 3.0tiA ± 0.9 MeV [127] ) although it could decay 
strongly into a pair of D mesons. Furthermore, the ir + ir~ invariant mass distribution 
first [123] and a detailed angular analysis next [122] showed that the dominant decay 
is X — > J/tpp, which would be isospin violating if the X(3872) were a conventional 
charmonium state. 
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Figure 7.: Invariant mass spectrum of the J/ipTt + ir system in the observation paper [123] 
(left) and of the D* D system in Ref. [127] 



The above-mentioned angular analysis from the CDF experiment [122] was able to 
discriminate among the possible J PC assignments, excluding any other possibility than 
J PC = 1 ++ and 2 h . In the meanwhile there came evidence for the decay X —> 
J/4>1 [128] and not of the decay X — > Xcil [123], thus confirming positive intrinsic 
charge conjugation. 

For several years the most favored option has been, since it has L = 0, J PC = 1 ++ . 
While we were finalizing this review, BaBar published a study of X — > J/ipuj where 
this option is disfavored with respect to the J PC = 2 h one, since it has only 7.1% 
probability to match data [129]. Implications of this result are still under study and 
cannot be reviewed here. 

4'1.2. Mass and width. The most recent advances on this topic concern the measure- 
ment of the mass of the X(3872) and the discussion on whether there is more than one 
state with similar mass (as predicted by tetraquark models) or not. The two options 
being investigated are that either the neutral and charged B mesons decay to different 
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Figure 8.: Measured mass of the V(3872) particle. The different production modes (B° — > 
XKg and B~ — > XK~ ) and the different decay modes (X — > J/tpirir and X —} D*°D°) 
are separated. 



linear combinations of the two possible X states or that the two states decay into different 
final states (in particular J/ipniT and D*°D°). 

The first possibility (different B a and B ± decays) has been investigated by CDF [137], 
where the J '/ipmr spectrum has been fitted searching for evidence of multiple structures. 
The negative result of such a search has allowed to establish that the eventual two 
states would have a mass difference smaller than 3.2 MeV at 90% C.L.. BaBar and 
Belle have instead measured the masses of the states observed in B° and B ± decays 
separately [130, 131], arriving at similar conclusions, AM = 1.2 ± 0.8 MeV. 

Both scenarios of multiple X states have instead been investigated by several mea- 
surements of mass of the X state performed by BaBar [130, 127] and Belle [131, 138] 
. The summary of all available mass measurements is shown in Fig. 7 where the mea- 
surements are separated by production and decay channel. The current word average is 
M = 3872.2 ± 0.4 MeV. There is an indication that the particle decaying into J/i/jmr is 
different from the one decaying into D*°D°, their masses differ by about 3.5 standard 
deviations. 

In addition, BaBar also measured the V(3872) width, T =(3.0±^ ± 0.9) MeV [127] 
, result which will be interpreted in the next section. 

4T.3. Production and decay. The X meson has been searched in several of its possible 
decay channels (/), by looking for B —> XK, X — >• /. The searched decay modes and the 
measured product branching fractions are listed in Tab XII. Representative measured 
spectra where the signal is present are summarized in Figs. 7 and 9. 

Several aspects of these measurements can be stressed: 

• albeit with relatively low statistics, the decay branching fraction into XK of the 
charged and neutral B decays are consistent: BF(B —> XK°)/BF(B —> XK ± ) = 
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Figure 9.: Representative measured invariant mass spectra of the most relevant de- 
cay modes where the X mesons has been searched. From left to right, top to bottom: 
[133], V(2S) 7 [133], y c (1P) 7 [123], J/^ [136] 



0.63 ± 0.16 for X ->■ J/^tttt and BF(B XK°)/BF{B XK±) = 1.5 ± 0.4 for 
X -» D*°D° 

• the DD* decay mode is the favourite by almost an order of magnitude 

• the ip(2S)~f branching fraction is comparable with the J/ip"f one. 

• the J/tprj decay has been searched for although it would violate charge conjugation. 
It was actually one of the evidences that lead to the determination of C. 

The measured product branching fractions can be translated into absolute branching 
fractions of the X particle by exploiting the upper limit on B — > XK measured by BaBar 
from the spectrum of the kaons recoiling against fully reconstructed B mesons [139], 
BF{B ± -» K ± X{'dS72)) < 3.2 x 10~ 4 at 90% CL. Combining the likelihood from the 
measurements of the product branching fractions in the observed channels, the B — > XK 
upper limit and the measured X width [127, 123] , with a bayesian procedure we extracted 
the likelihood for the absolute X BF and the widths in each of the decay modes. Then, 
we used the probability distributions extracted from this procedure to set limits on the 
not observed channels. The full shape of the experimental likelihoods was used whenever 
available, while gaussian errors and poissonian counting distributions have been assumed 
elsewhere. The 68% confidence intervals (defined in such a way that the absolute value of 
the PDF is the same at the upper and lower bound, unless one of them is at the boundary 
of the physical range) are summarized in Tab. XII for each of the decay modes. Fig. 10 
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shows the likelihoods for the X branching fraction in J/tpTrir and D*°D°, the total width 
and the partial width T(X -> D*°D°). 
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Figure 10.: Likelihood function of the X branching fraction in J/ipirir and D*°D° , the 
total width and the partial width T(X — » D*°D°). See text for a description of the 
combination method. The dark (light) filled area corresponds to the 68% (90%) C.L. 
region. 



Finally, information can be obtained also on the production mechanism. The likeli- 
hood combination described here allowed also to extract the B — > XK branching frac- 
tion. The result (see Fig. 11) shows that the X meson is less copiously produced than 
all other charmonia, for which the corresponding branching fration is at least 5 x 10~ 3 . 
As far as other production mechanisms are concerned, both the B — > XK* decays and 
the 77 production have been investigated. The measured BF(B° X(3872)K* (892)°) 
shows that the X meson is produced in association with a K* less favorably than in 
association with the K, contrarily to the other charmonia. The only search in 77 pro- 
duction was performed by CLEO [140]. It did not return a signal and a limit was set: 
(2J + l)r 77 fl(X- > 7T + 7T - J/ip) < 12.9 eV. 

4'2. The 3940 family. - Three different states have been observed in the past years 
by the Belle collaboration with masses close to 3940Mev/c 2 . Their measured masses and 
widths are summarized in Tab. XIII and Fig. 12. 

The first one, named X, was observed in continuum events (i.e. not in Y(4S) decays) 
produced in pair with a J/ip meson [144]. Subsequently its decay into DD* has been 
ascertained by means of a partial reconstruction technique in the same production mecha- 
nism [141]. The production mechanism constrains it to have positive charge conjugation. 
Furthermore the fact that in this production mechanism only J = states are observed 
makes it likely to be either J PC = 0++ or 0~ + . Ref [144] was also able to measure the 
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BF(B ->XK) 

Figure 11.: Likelihood function of the B — > X branching fraction, compared with the 
measured branching fraction for the corresponding decays with regular charmonia. See 
text for a description of the combination method. The red (yellow) interval corresponds 
to the 68% (90%) C.L. region. 



absolute production rate of this state and searched for X — > J/ipcu without evidence of 
signal. It eventually concluded that BF(X -> DD*) > 41% and BF(X -» J/ipco) < 26% 
90% C.L.. 

A second one, named Y, observed in B decays and decaying into J/ipu; [53], is sig- 
nificantly larger, and only its charge conjugation (C=+) is known, while J = 0,1,2 
and all parities are possible. The BaBar collaboration confirmed the F(3940) — > J/iJjuj 
decay [54], but measuring a lower mass and a width, albeit marginally consistent (see 
Tab XIII) . While the experimental differences among the two measurements are still un- 
der investigation, a new study has observed a state, named F(3915) which is consistent 
with the parameters as measured by BaBar, but produced in 77 collisions and decay- 
ing into J/ijju [142] (see Fig. 13). Under the hypothesis that the F(3940) and F(3915) 
states coincide, this is the first case of two production mechanisms being observed. The 
y(3940) — » DD* decays have been searched in Ref. [138]. In the absence of a signal a 
limit on B_F(F(3940) -> J /i/)U))/BF(Y (3940) DD*) > 0.71 90% C.L. is set. This 
result supports the hypothesis that this is a different state from the X(3940), where 
BF(X(3M0) -> J/iJ;uj)/BF(X(3M0) -> DD*) < 0.58 90% C.L.. 

The third observed state, named Z, is produced in two-photon reactions and decays 
into D-pairs [143]. This state has been confirmed in Ref. [145], where the J PC = 2 ++ 
has been ascertained. 

Because of their quantum number assignments and their masses these states are good 
candidates for the radial excitation of the x mesons, in particular the Z(3940) meson 
could be identified with the Xc2^P) and the y(3940) with the x c i(2P). The unclear 
points are the identification of the X(3940) state and the explanation of why the F(3940) 
state does not decay preferentially in D mesons. 
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Figure 12.: Measured masses and widths of the "3940 family" of states. 



4'3. Other C — + states. - Another set of states has been observed with C — + and 
masses close to 4150 MeV, creating another interesting cluster of states. 

The first of such states was observed by Belle produced in continuum events and in pair 
production with a J/ip meson and decaying into D* + D*~ [141] (see Fig. 14. Its measured 
mass and width are M = Al56±f (stat.)±15(sys.) MeV and T — 139±l\ 1 (stat.)±21(sys.) 
MeV. The production mechanism favors J = and it enforces C = +, making it a good 
candidate for rj c (3S), which would have J PC = CT + . 

Next comes the report from CDF of an enhancement close to the threshold in the 
J/ip(f> invariant mass in B — > J/iptfiK decays [59]. The fitted mass of the resonant state, 
Y(4140), is M = (4143.0 ± 2.9 ± 1.2) MeV with T = (11. 7^ 8 3 ± 3.7) MeV. While the 
natural (L — 0) quantum numbers assignment would be J PC = ++ , J PC = 1 ^ is 
also allowed, leaving open the possibility of this state being the lowest lying hybrid. This 
hypothesis is also suggested by the closeness of the measured mass with the predictions of 
lattice calculations for the hybrid ground state (see for instance Ref. [146]). The search 
for such a state in this production mechanism at _B-Factories does not have enough 
statistics to be conclusive. Belle has therefore searched this state in 77 production which 
is expected to be copious in certain models [147]. No evidence was found and a limit 
r 77 x BF((}>J/ip) < 41 (6) eV for J p = 0+ (2+) were set 90% C.L. for the Y(4140). 

When performing this search a 3.2tr evidence was found with M = 4350.6l5;'i i (stat) ± 
0.7(syst)) MeV and T = (13±g 8 (stat) ± 4(syst)) MeV. Yet another state to place in the 
spectrum, the V(4350), with C=+, and close to one of the J PC = 1 states reported 
in Sec. 4'4. 

4'4. The 1 family. - The easiest way to assign a value for J PC to a particle is 
to observe its production via e + e~ annihilation, where the quantum numbers must be 
the same as the the photon: J PC = 1 . B factories can investigate a large range of 
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Figure 13.: Invariant mass distributions of the most significant observations of the "3940 
family" states: X(3940) -> DD* , F(3940) J/^w, F(3915) J/V>w and Z(3940) -> 



masses for such particles by looking for events where the initial state radiation brings 
the e + e~ center-of-mass energy down to the particle's mass (the so-called 'ISR' events). 
Alternatively, dedicated e + e~ machines, like CESR and BEP scan directly the center- 
of-mass energies of interest. 

The observation of new states in these processes started with the discovery of the 
F(4260) — > J/t/jir + ir~ by BaBar [70], promptly confirmed by CLEO-c both in the same 
production process [148] and in direct production [149]. The latter paper also reported 
evidence for F(4260) -> J/ijm ?: and some events of F(4260) -> J/ipK+R-. 

While investigating whether the F(4260) decayed to ■0(2iS , )7r + 7r~ BaBar found that 
such decay did not exist but discovered a new 1 state, the y(4350) [79]. While the 
absence of F(4260) — > ip(2S)ir + ir~ decays could be explained if the pion pair in the 
J/tp7r + ir~ decay were produced with an intermediate state that is too massive to be 
produced with a ip(2S) (e.g. an /°), the absence of F(4350) — > J/ipTr + ir~ is still to 
be understood, more statistics might be needed in case the y(4260) decay hides the 
y(4350). 

Next, Belle has published the confirmation of all these 1 states [78, 81] and at the 
same time has unveiled a new state that was not visible in BaBar data due to the limited 
statistics: the y(4660). Fig. 15 shows the published invariant mass spectra for both the 
J /ipir + ii~ and the ■0(2<S')7r + 7r - decays. 

Critical information for the unravelling of the puzzle is whether the pion pair comes 
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Figure 14.: Invariant mass distributions of the most significant observations of the states 
with C = + and mass above 4 GeV: X(4160) -> D*+D*- ; X(4140) -> J/^w, and 
X(4350) -> ZXD. 



from a resonant state. Fig. 16 shows the di-pion invariant mass spectra published by Belle 
for all the regions where new resonances have been observed. Although the subtraction of 
the continuum is missing, there is some indication that only the Y (4660) has a well defined 




Figure 15.: J/tJjtt + tt (left) and %p(2S)n + TT (right) invariant mass in ISR production. 
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Figure 16.: Di-pion invariant mass distribution in y(4260) — > J/iJm + Tt (left), 
F(4350) -> iP(2S)tt + tt~ (center), and F(4660) -> ^j(2S)-k + tt~ fngfttj decays. 



intermediate state (most likely an /o), while others have a more complex structure. 

A discriminant measurement between Charmonium states and new aggregation forms 
is the relative decay rate between these decays into Charmonium and the decays into 
two charm mesons. Searches have therefore been carried out for Y -> £>(*)£>(*) de- 
cays [150, 151, 152] without any evidence for a signal. The most stringent limit is [152] 
BF(Y (4260) -> DD)/BF(Y (4260) -> J/^tt+tt-) < 1.0® 90% confidence level. 
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Figure 17.: The A C A C distribution in ISR (left) and the J/^jlu distribution in 77 events 
as measured by Belle. The superimposed lines are the result of the fits to the data whose 
results are reported in the text. 



A distinctive signature of tetraquarks would instead be the observation of decays of 
these mesons into two baryons. This has lead to the search by Belle of resonant structures 
decaying into A C A C associated to ISR [80] . A tantalizing structure appears near threshold 
(see Fig. 17, which fitted with a Breit-Wigner returns M = (4634 j^tg) MeV and T to t = 
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(92^24^12) MeV, properties which are close to the one of the y(4660). A consistent 
analysis of the two measured spectra has actually concluded that the two structures 
come from the same state with a large preference for the baryonic one (BF(Y(4660) — > 
A C A C )/BF(Y(4660) -> xp(2S)%-!r) = (25 ± 7) [153]). 
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4'5. Charged States. - The real turning point in the quest for states beyond the 
Charmonium was the observation of charged states decaying into charmonium. There is 
in fact no way to explain such observation without having at least four bound quarks 
(ccud). There is currently evidence of three charged states, all seen exclusively by the 
Belle Collaboration: the Z(4430) state decaying into , 4)(2S)tt ± [154], and the Z\ and Zi 
states decaying into Xci(2S)n ± [57]. 

Unfortunately these states have been observed in B decays in association with a 
charged kaon, i.e. in three body B — > X cc ttK decays, where X cc — ip(2S) or Xd- Three 
body decays suffer from interferences between strong amplitudes mediated by different 
resonances. In these particular cases the ttK system presents several known resonances 
that could cause significant effects of reflection. Namely, the decays B X CC K* (892), 
B —> X CC K* (1410), and in particular their interference constitute irreducible sources of 
background which are difficult to estimate. The original observation of the Z(4430) [154] 
has been therefore object of detailed scrutiny by the BaBar collaboration who performed 
the search for the same final state by studying in detail the efficiency corrections and the 
shape of the background taking the latter from the data as much as possible [155]. The 
search resulted into hints of a structure close to Belle's observation, but after accurate 
estimate of the background results into an exclusion on the product of branching fractions 
BF(B Z (4330) K + )BY(Z (4330) ip(2S)n) < 3.1 ■ 10~ 5 @95% C.L. to be compared 
with BF(B -> Z (4330) K+)BF(Z (4330) -> ^(2S)tt) = (4.1^;°) • 10~ 5 as reported in 
Ref. [154]. Above all the paper from BaBar has raised the attention to the problems 
inherent to the analysis of the Dalitz plot, leading to a reanalysis from Belle [156]. 

The results of the latter analysis, that confirmed the original ones although the errors 
on the parameters increase significantly, and of the B — > Xci^K decays (not yet analyzed 
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with the full three body approach) are shown in Fig. 18. The fits to the ip(2S)ir and 
Xci7r invariant mass distributions return masses M = 4443^8 MeV, M = 4051 ± 14rt|5 
MeV, and M = 4248±^|±|| a MeV, and widths V = 109l^ 3 MeV, V = 82+^+^ MeV, 
and T = 177±|f±|i 6 MeV for the ^(4430), Z x , and Z 2 states respectively. 




Figure 19.: Two more representative spectra that have not been analyzed and that will be 
discussed in the concluding section of the paper: pp invariant mass spectrum from B — > 
ppK decays as in Ref. [157] and D*°D*~ invariant mass spectrum from B — > D D*~K 
decays as in Ref. [158] 



4'6. Interpretations . - In this section we match the information of the observed states, 
summarized in the previous sections, with the possible theory interpretations detailed 
in Sec. 2, we detail the arguments in favour or against the different interpretations and, 
when possibile, mention the measurements that would help resolve the ambiguities. 

4'6.1. X(3872). V(3872) represents a very peculiar case because it lies very close to 
the DD* threshold, its binding energy with respect to this threshold being compatible 
with zero. This characteristic lead many authors to identify the X with a DD* molecule 
at threshold. Since nix < rn D o + m D o* , at least in the J/i/>7r + 7r~ channel, the state 
could be considered as a proper bound state, i.e. with negative binding energy (Eg = 
m x — mrjo — m^o* < 0). 

This interpretation, as the others we will review in the following, stand on the hy- 
pothesis that the X has J PC = 1 ++ . This prejudice came from the fact that even if 
CDF stated in [122] it could not distinguish between 1 ++ and 2 ^ , the Belle analysis 
[159] seemed to prefer the 1 ++ assignment. 

The existence of such a bound state has been investigated using the potential models 
discussed in Sec. 2'2. A lot of work has been spent since the first observation of the X 
to obtain solid and quantitative results from those potential models, but yet there is not 
a final answer in the literature to the question "Does a 1 ++ DD* bound state exist at 
threshold?" . 

Tornqvist was the first to state in [29] and later to confirm in [160], that a DD* bound 
state for J PC = 1 ++ and / = is likely to exist. From the point of view of one pion 
exchange the problem is tightly related to the deuteron case. In the deuteron case the 
tensor term of the potential in Eq. (46) favors the binding allowing for S to D wave tran- 
sitions. Indeed the scalar part of the potential provides only about one third of the total 
binding energy. It is important to notice that one cannot predict exactly the deuteron 
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binding energy without a detailed knowledge of the very short range interactions, since 
its value cannot be deduced from the other scales of the problem. Nevertheless once the 
pion-nucleon coupling constant is known and a regularization scheme for the potential is 
given, one can use the measured binding energy to deduce the value of the cutoff A and 
thus the long distance part of the wave function. The same reasoning applies to the X 
case. The binding channels for a DD* state are J PC = h , 1 ++ with 1 = 0. Assuming 
some definite value for the binding energy, one is able to compute the wave function at 
long distance whose components are D°D°* and D + D~* both in 5-wave and f-wave. 

A slightly different result is found by Swanson in [25], where the combined effects 
of short-range one-gluon-exchange potential and long-range one-pion-exchange potential 
are required to obtain the binding. The possibility to include other channels, such as 
J/i/jp, J/tpu> and D + D~* is investigated. 

On the other hand, a completely different result is obtained by Suzuki in [161]. The 
effective potential between D and D* due to one pion exchange is computed from the 
Born scattering amplitude: A = (e'* ■q)(e-q)g 2 / (rn^ — q 2 — ie), where q is the it momentum 
and thus the difference between the D and D* momenta. Since the decay D°* — > D°tt° 
occurs with a tiny Q-value A = to^* — mjj— ~ 8 MeV the pion is produced practically 
at rest and thus in the non relativistic limit for the denominator of the pion propagator 
is: m 2 — q 2 ~ — (2m. K A — |q| 2 ). When computing the Fourier transform of the scattering 
amplitude one obtains a three-dimensional delta function plus correction of C((2m w A) 2 ), 
which can be neglected. Unlike the one-dimensional ^-function, the three-dimensional S- 
function in a Schroedinger equation does not admit bound states. The author concludes 
that the X cannot be regarded as an analog of the dcuteron and furthermore questions 
that the binding is recovered when the short-distance one-gluon-exchange potential is 
added. Actually for a spatially extended object, as the X is, it is difficult to believe 
that the binding is provided by short range forces. Furthermore the possibility that 
the binding is due to the coupling with other molecular channels, such as J/i/joj, seems 
very unlikely since the coupling to this channel would occur through D/D* exchange 
thus resulting in a Yukawa potential with very short range (^0.1 fm). On the other 
hand Suzuki considers more likely that the binding proceeds through the coupling to 
charmonium states, due to u,d quark exchange. However in this case the nature of the 
state would be essentially that of a charmonium and the dynamics completely different 
from the molecular one. For a recent review on the one pion exchange binding mechanism 
see [162]. 

Besides the molecular interpretation, alternative explanations have been proposed, 
such as a regular charmonium assignment [163, 164] and a tetraquark interpretation 
[3, 4]. We will review the characteristics of the different interpretations by topic. 

Quantum numbers. From the point of view of quantum numbers many considerations 
were made soon after the first observation of X in 2003 [165]. The absence of the decay 
in D°D° despite the ample phase space (decay momentum ~500 MeV) suggests 
unnatural spin parity J p = 0~ , 1 + , 2~, . . . or unnatural spin-charge conjugation parity 
J PC = H , 1 h , 2 H , . . .. The dipion quantum numbers help discriminating among the 
possible J PC of the X. A scalar dipion (J PC = ++ ) combined in S-w&ve with the J/ip 
gives a 1 state, in F-wave a 1 H state and in D-w&ve a 2 state, all with isospin 
J = 0, 2; a vector dipion (J PC = 1 ) in 5-wave with the J/ip gives a 1 ++ state, while 
in P-wave a 2 ^ state, all with isospin 1=1. Finally a tensor dipion (J PC = 2 ++ ) gives 
in 5-wave a 2 state with isospin I = 0,2. Both the J PC = 1 ++ and 2 ^ assignments 
favor the vector dipion configuration, which would also explain the predominance of the 
decay into J/ipp. 
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As stated in the beginning, the vicinity of the X mass to the DD* threshold has 
prompted considerations of the possibility that X could be a DD* molecule: 



DD*) ± \DD*) 

Potential models predict that 1 = states are favored with respect to / = 1 states, in par- 
ticular they predict h and 1 ++ bound states. This would imply an equal contribution 
of the charged and neutral mesons components: 



(61) 



\Z)0 D o*\ + \ D + D -*\ 
(62) \DD*) = J >±i > 

However since the binding energy is much smaller than the mass gap between neutral 
and charged pairs, one expects substantial isospin breaking in the wave function. Indeed, 
since the D+*D- threshold (-3879 MeV) lies -8 MeV above the D°*D° threshold 
(~3871 MeV) it should have a smaller weight than D°*D°, which means that there will 
be a strong 1 = 1 component in the state. Some new results on the isospin issue have been 
obtained in [166, 167, 168] studying the charged versus neutral D mesons components of 
the X wave-function. 

As for the charmonium option two assignments have been proposed [163, 164, 169]: 
ID and 2P respectively for a J p = 2" and a J p = 1+ state, since the parity and charge 
conjugation quantum numbers of a pure cc state are related to spin and orbital angular 
momentum as P = (— l) i+1 and C = ( — l) i+5 . A computation of the radiative decays 
and of the principal hadronic decay modes allows to exclude some of the states among 
ID and 2P, due to the extremely small width of the X. The states which remain left are 
1 3 -C>3,2, 1 1 P ) 2, 2 3 Pi and 2 1 P\. Since the X has even charge conjugation, there remain 
two possibilities: 1 1 Z?2 and 2 3 Pi, respectively for 2 h and 1 ++ J PC quantum numbers. 

The main limit of the charmonium assignment, which seemed to be totally ruled out, 
is the isospin quantum number. The observation of the decays X — > J / ipir + ir~ ir° and 
X — > J/ipTr + Tr~ with almost equal branching ratios indicates that X contains an equal 
amount of I = 1 and 1 = components. A standard cc state cannot account for this, since 
it is a pure isoscalar state. Colangelo and collaborators argued that the isospin violation 
is not so severe [170] due to the different phase space volumes available respectively for 
J/ipTT + TT~ and J/ipTT + ir~Tr° final states. This is due to the different total decay widths 
of co and p. Taking this into account the 1 = 1 amplitude is approximately five times 
smaller than the 1 = one. This smaller isospin violating effect could be explained with 
the difference in mass between the charged and neutral D-mesons as intermediate states 
in the X decays. 

Finally the tetraquark. In this approach the 1 ++ diquark-antidiquark bound state, 
with symmetric spin contribution: 



(63) X q = [cq]s=i[cq]s=o + [cq]s=o[cq]s=i 

is a candidate to explain X, since it is expected to be narrow like all diquark-antidiquark 
systems below the baryon-antibaryon threshold. Furthermore unnatural spin parity for- 
bids the decay to D°D°, while the decay to J/ip plus vector meson is allowed with 
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conservation of the spin of the heavy quark pair. The isospin quantum number is re- 
lated to the finer structure of the X q state. The two flavor eigenstates X u and Xd mix 
through self energy diagrams, which annihilate a uu pair and convert it into a dd pair 
through intermediate gluons. In the basis {X u , Xd} the annihilation diagrams contribute 
equally to all the entries of the mass matrix, while the contribution of the quark masses 
is diagonal. The resulting 2x2 mixing matrix is: 

<«> ("Y' J +s ) 

At the scale determined by the cc pair the annihilation term 5 is expected to be small 
and thus the mass eigenstates should coincide with flavor eigenstates to a rather good 
extent. One can put: 

Xi ow = cos 9X U + sin 9X d 

(65) 

Xhigh = - sin 9X U + cos 9X d 

where, in the limit 5 << m u ^d, cos# ~1 and sin# ~ 9, giving an almost maximal isospin 
breaking. The mass difference between the two states is: M(Xhi g h) — M(Xi ow ) = (m<j — 
m u )/ cos29. Indeed, as already mentioned, the possibility that there were two different 
states A(3872) and X(3875) is not excluded from experimental data. An alternative 
mechanism to explain the isospin violation in the tetraquark picture is the possibility of 
au-p" mixing, as proposed in [171]. 

Mass and width. Beside the vicinity of the X mass to the D°D°* threshold, one of 
the characteristics of the X is its extremely narrow width, that needs to be accounted 
for. 

In the molecular picture the mass of the X is directly related to the D°D°* threshold. 
For what concern the width of the state the main arguments follow from the low- energy 
universality [172] of such a loosely bound molecule. Indeed an 5- wave D°D°* molecule has 
an extremely small binding energy E B , which results in an unnaturally large scattering 
length a = ^j2\iE B . Therefore E B =0.25 MeV gives a ~ 8 fm, which is larger than the 
typical range of strong interactions. Thanks to this feature the molecule has properties 
that depends on a, i.e. universal properties, insensitive to the details of the interaction 
between mesons. Low-energy universality implies also that the asymptotic form of the 
D°D°* bound state wave function is known and can be expressed in terms of a: ip{r) = 
(l/v / 2~7ra)e~ r / a /r. In addition to the D°D°* component there will be contributions from 
the other 1 ++ hadronic states H with nearby threshold. In general one can write: 

(66) \X) = Z%, (\DD*) + \DD*)) + £ Z^H) 

H 

Low-energy universality implies that as the scattering length a increases, the proba- 
bilities for states other than D°D°* decrease as Zh ~ 1/a and in the limit a — > oo the 
state becomes a pure D°D°* molecule. It is important here to notice that the X rep- 
resents a really particular case of bound state, since its binding energy is much smaller, 
i.e. compatible with zero, than the natural energy scale m\j ' [i ~20 MeV. 
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In this framework one can compute the partial widths in the different decay channels. 
As for the constituent decay channels one has: 

T(X -> D°D°~f) = Znn^TiD * -> £° 7 ) 

where CV and C 7 take into account the interference effect between the D°* and D°* 
decays. This interference effect was computed for the first time by Voloshin in [173] 
and is able to account for an enhancement factor of the width of 2 at maximum. As 
for the short-distance decay modes (J/ipp, J/tpco, J/tpl, V'(2'S , )7), since Zh ~ 1/a, the 
main contribution will come from the D°D°* component of the wave function. One can 
estimate that these decay rates scale like 1/a as a — > oo. In this way the molecular 
picture accommodates the narrowness of the X. We will discuss some problems of this 
picture in the last paragraphs of this section. 

As for the charmonium assignment potential models predict masses which are smaller 
than 3872 MeV for the ID case, while larger than 3872 MeV for the 2P one [163]. The 
l 1 !^ state should have a quite large radiative branching fraction to XcT- B(\ l D2 — > 
Xcl) ~ 0.5, which instead has not been observed until now. On the other hand the 
2 3 Px charmonium is predicted to have branching fractions of a few percents to J/ipj 
and ip(2S)^, which seems in agreement with Tab. XII. However, in [164] a study of 
the effect of the open charm thresholds on the charmonium properties gives very large 
partial width into D°D°* for 2 3 Pi which are not compatible with the narrowness of the 
X. The \ x Di state on the contrary is predicted to have a large branching fraction to 
h c j, which has not yet been looked for. The experimental search for this decay mode 
would be thus of crucial interest. The interplay between the open-charm mesons and the 
charmonium states has been indicated as a possible solution for the identification of the 
X. Many authors computed the effect of virtual or real D- mesons loops in shifting the 
masses of ordinary cc states, for a recent review see [174]. Of particular interest for the 
X case is a study by Kalashnikova [175], in which a bound state with J PC = 1 ++ nearly 
at threshold is obtained from the interaction of the 2 3 P X charmonium level with DD* 
mesons. For an update on this subject see [176]. 

In the tetraquark picture, as explained in [3, 4], the mass of the 1 ++ candidate state 
can be written in terms of the chromo-magnetic coupling and the [eg] diquark mass. 
Taking for the couplings the values obtained in [3] and assuming as an input the mass of 
the X one can deduce the mass of the [eg] diquark: m[ cg ] =1933 MeV. From this value 
the entire spectrum of [eg] [cq] can be computed [3] . The finer structure of the X q states 
has been already discussed and leads to the prediction of two different states, which 
can be identified, as a first approximation, with the two flavor eigenstates X u and X4. 
The observation in 2006 of a state decaying to D a D°ir with mass 3875 MeV favored the 
assignment: X u = X (3875), decaying mainly into J/t(nr + ir~ and Xd = X(3872) decaying 
into D°D°tt. The mass ordering of these two neutral states seems to be reversed, since 
the u quark is lighter than the d quark and thus one would expect X u to be lighter 
than Xd- However the quarks which form the diquarks in the X u have the same electric 
charge and thus a consistent consideration of the electrostatic energy can perhaps change 
the order of the masses. Besides these two neutral states, two charged states arise as a 
natural prediction of the tetraquark picture X + = [cu][cd] and X~ = [cd][cu]. The lack 
of any observation of these two states constitutes the main drawback to the tetraquark 
assignment. 
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Production and decay. The production of a D°D°* molecule necessarily proceeds 
through the production of a D meson pair. This production mechanism has been exten- 
sively studied in B decays by Braaten and collaborators in [177]. For a bound state at 
threshold standard quantum mechanics predicts that the wave function in momentum 
space has a Lorentzian shape. From this information one can extract the elastic scat- 
tering amplitude A(D°D°* -> D°D *) and the transition amplitude A(D°D°* X) as 
the residue at the pole of the elastic scattering amplitude. 

This allows to compute the partial width of B ± — > XK ± and the differential dis- 
tribution aT(B ± — > D°D°* K + ) /dMjyjj* as functions of model dependent parameters, 
describing the formation process of the X as a coalescence process between D° and D°* 
mesons. The distribution oT(B ± D°D°*K+)/dM DD » shows a peak more and more 
evident near the threshold as the binding energy is decreased. In this computation the 
authors make use of an effective held theory description of the interaction between the 
two £>-mesons. The effective field theory method has been first introduced in [173] and 
then used in [172, 178] to compute some properties of the X meson, relying on the fact 
that as far as the exchanged momentum between the Z?'s is small, they can be described 
with elementary fields. In this approach a coupling to the regular charmonium state has 
also been implemented [172]. 

In [179] the exclusive production of X in B decays is studied. The most relevant 
result of this paper is the prediction of a deviation from the standard charmonium case 
of the ratio K = B(B° XK°)/B(B+ XK+) ~ 0.1 (in the charmonium case 
1Z <~ t(B°)/t(B + ) <~ 0.9). The BaBar collaboration has measured [130] the two products 
of branching ratios: B(B+ -> K+ X) x B(X -> J/ipn+n-) and B(B° -> K°X) x B(X -> 
J/i(jTr + iT~). These measures give 1Z — 0.6 ± 0.2. Instead for a tetraquark state 1Z = 
0.5 ±0.3 [4]. 

In the molecular picture the decay into J/t/j + light hadrons has been accounted for by 
including the possibility for inelastic channel in the D°D°* scattering introducing a com- 
plex scattering length [180]. The authors find factorization formula for the decay width 
separating the long distance term (the transition amplitude X — > D°D *) characterized 
by a typical distance of a ~ 8 fm, from the short distance one, acting on distances of 
about 1 / nij, . The first term can be computed as before as the residue at the pole of the 
elastic scattering amplitude and does not contains any unknown parameter, while the 
second one is related to the coupling of the X to the hadronic channel considered. This 
is estimated from the wave-function component associated to each hadronic channel and 
can be borrowed from the potential model results. 

Finally studies of the long-distance D°D°ir° decay mode has been carried out to 
obtain the relative partial width and the modification of the line shapes of the resonance 
due to the vicinity to the threshold. The study of the line shapes in this channel has 
been useful to distinguish between a bound state below threshold and a virtual state 
above threshold. Indeed the temporary agreement between Belle and BaBar data on the 
evidence for an higher mass value of the peak in the D°D°* invariant mass spectrum, 
lead many authors to solve this problem identifying the X with a virtual state. In [181] 
it was performed an analysis on the Belle data [132]. This analysis concluded indeed 
that the resonance observed at 3875 MeV in the D°D°* channel could be related to the 
one observed at 3872 MeV in the J/iprnr channel only if the X is taken to be a virtual 
state. The consequences of this conclusion for the isospin breaking have been considered 
in [182]. In a subsequent analysis [138], with higher statistics, the Belle collaboration 
claimed that the resonance peak in the D°D°* channel appears at <~ 3872 MeV. In 
[183] the authors perform a re-analysis of the new data from Belle [138] together with 
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the ones from BaBar [127] on the X decays to D a D°ir. There are two main novelties in 
this work: (1) they include a finite width for the D°* constituent; (2) they consider the 
fact that the invariant mass of a D°tt° pair coming from the decay of a D°* in the X 
is not equal to the mass of a D°* meson since the energy of a D°* meson bound inside 
the molecule differs from its mass by about the X binding energy. They derive the line 
shapes of the resonance in the two final states and fit the two data sets separately. They 
obtain the position of the peak, its width and some ratio of branching fractions. The 
fit to the D°D°ir decay mode returns a real scattering length, which means that the 
transition amplitude to the inelastic channels (in this case J / ipiT + ir~ ) is null. From the 
same fit they obtain that B(X -> D°D°ir°)/B(X -> J/^tt+tt-) = 0.004, which instead 
implies a strong predominance of the inelastic decay channel ( J/'i/>7r + 7r~) over the elastic 
one (D D tt°). The authors adduce this inconsistency to an improper treatment of the 
experimental resolution. 

The same analysis has been repeated and upgraded in [184, 185]. In particular [184] 
included the ip 7 decay mode. The results are consistent with X being a true bound 
state. The line shapes have been studied also in [186], where a general discussion about 
near to thresholds pole in the S-matrix is carried on, and in [187], where the authors state 
that the X resonance is the result of a fine-tuning between the regular 2 3 Pi charmonium 
state and the D°D°* threshold. 

Other decay modes have been considered in the hypothesis of X being a loosely bound 
D°D°* molecule. In particular Voloshin and Dubynskiy [188, 189] studied the transition 
X ~~ XcjTT , considering the cc pair inside the X to be in a spin-triplet configuration, 
since it is favored by even charge conjugation. In [190] the X — » DD-f decay mode is 
taken into account, since the photon spectrum would provide information on both the 
long and short distance component of the molecular state. 

Concerning the production mechanism in pp collisions, as we have already stated in 
Sec.2'2, the production of a loosely-bound S-wave molecule must proceeds through the 
production of D-mesons pairs with suitably small relative three-momentum to account 
for the tiny binding energy. While in the B-mesons decays the Z?-meson pair is produced 
practically at rest (Mb — Mk — Mn — Mo, ^900 McV), it seems odd that a such a 
loosely bound molecule could be produced promptly in a high energy hadron collision 
environment. Indeed one can learn some general features of the X as a molecule without 
any knowledge on the interaction details. The scattering length can be interpreted as 
the spatial extension of the hadron itself and is estimated to be <~ 8fm, much more than 
that of an ordinary hadron. Using a minimal uncertainty principle one is able to extract 
the relative momentum spread of the pair of mesons bound inside the X: Ak ~ 15 MeV. 
The central value for the relative momentum is of the order of the decay momentum of 
the D-mesons: k ~ 30 MeV. 

Starting from these considerations [191] computed the maximum theoretical value 
for the prompt production cross section of the X in pp collisions at the Tevatron. The 
authors make use of a Schwartz inequality, which allows to write the theoretical prompt 
production cross section as an integral over the differential cross section for the production 
of a pair of D°D°* mesons with fixed relative three-momentum. The integral extends up 
to ko + Ak ^50 MeV. The differential distribution for the cross section can be computed 
using standard MC event generators such as Herwig and Pythia. The results are shown 
in Fig. 20. The maximum integrated cross section amounts to erj^ ax — 0.085 nb for both 
the MC event generators. 

This estimate has been compared with a lower limit for the experimental cross sec- 
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Figure 20.: Left: the integrated cross section obtained with Herwig as a function of the 
relative momentum of the mesons in the D°D°* molecule. This plot is obtained after 
the generation of 55 x 10 9 events with parton cuts p^ art > 2 GeV and yP art > 6. The 
cuts on the final D mesons are such that the molecule produced has a p± > 5 GeV and 
y < 0.6. Right: same plot but using Pythia. 



tion obtained by the CDF collaboration properly taking into account the rapidity and 
transverse momentum cuts: the CDF collaboration recently performed an analysis [122, 
192, 193, 194] to distinguish the fraction of X and ip(2S) produced promptly from that 
originated from B-decays and found that only the (16 ± 4.9(stat.) ± 2.0(syst.))% of the 
X produced in pp collisions comes from B decays. Using the corresponding tJj(2S) result 
one can obtain an estimate of the product cr^ mpt (pp — > X + al\) x B(X — > J/?/>7r + 7r~) = 
(3.1 ± 0.7)nb. From the results on B(X — > J/ipTr + ir~), shown in Sec. 4T one can con- 
clude that a exp (pp — > all) < 33nb. Thus the upper bound on the theoretical cross section 
is ^300 times smaller than the lower limit on the experimental cross section, putting in 
trouble the molecular interpretation of the X. 

Nevertheless, in a subsequent paper by Braaten and Artoisenet [195] the effect of final 
state interactions (FSI) between the D° and D°* mesons has been taken into account 
to reconcile the experimental and theoretical estimates. According to the authors of 
[195] FSI should be able to rescatter higher relative momentum pairs to lower relative 
momentum pairs, thus allowing to consider k max ~ 2m w and thus to integrate a larger 
cross-section. Furthermore they introduce an enhancement factor in the cross-section 
~ l/(fc 2 + a 2 ), as stated in the Migdal- Watson theorem. In this way they succeed in 
reconciling the theoretical and experimental values. 

This result is strictly related to the use of the Migdal- Watson theorem, which formal- 
izes the FSI mechanism. There are two main conditions that have to be met in order 
to safely apply the theorem. First, the relative momentum between the two rescattering 
particles has to be smaller than the inverse of the range of strong interactions, namely 
200 MeV. Second, it is necessary that no more than two hadrons match the first con- 
dition. Reference [196] showed that the latter is not satisfied. Furthermore going up to 
relative momentum of ~ Im^ one is no allowed to ignore higher partial wave scattering 
beside the S-wave one. However if FSI works in the direction of an enhancement [196] 
predicts that if any X s = D S D* exists it should have an observable cross section at the 
Tevatron, namely cr^ ax (pp — > X s -fall) ~ (2± l)nb. It could decay to J/ipKK. This can 
be used as a test table for our understanding of FSI and of meson-meson interaction at 
the same time. 
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Finally another possible production mechanism for the molecular X has been proposed 
in [197, 198] to be e+e" ->• jX. 

As for the charmonium assignment one could refer to [199] for a complete review of 
charmonium production in high energy collisions. Concerning the J = 2 ^ possibility 
a study of gluon fragmentation to _D-wave charmonia is contained in [200]. 

In these final paragraphs we discuss some new considerations about the D a D°Tr a 
partial decay width and the radiative decay modes. We think that these two arguments 
may offer an additional way to test the supposed molecular nature of the X. 

From the analysis described in Sec. 4*1, which summarizes all the experimental 
informations available up to now, we are able to extract a probability density function 
(PDF) of the partial width of the X -> D°D°* -> D°D°tt (see Fig. 10) and consequently 
estimate the probability T(D°D *) <100 KeV to be - 0.9% 



o. 



This result can put in some trouble the molecular interpretation of the X. Indeed, 
due its small binding energy one would expect the D°* and D°* mesons bound inside the 
X to decay as if they were free, that is with the decay widths measured for free mesons. 
In the following we will refer to the decay D°* — > D°ir, but the same considerations are 
valid for the charge conjugate mode. We will discuss in the end the effect of interference 
between the charge conjugate modes. 

The measured partial width for D°* -> D°ir is T(D°* -> D°tt) = 70 KeV. Neverthe- 
less when considering this decay mode for a meson bound with a D° meson to form the 
X we find ourselves in a rather peculiar situation: the process D°* — > D°ir is the same 
process which is responsible for the binding of the molecule. It has been claimed by many 
authors that the one pion exchange potential is sufficiently attractive in this channel to 
admit a bound state exactly at threshold. Thus it can occur that a would-be-real-pion is 
reabsorbed by the other meson keeping the molecule bound and making it living longer 
than its constituents. The relevant diagrams are reported in Fig. 21. 
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Figure 21.: Feynman diagrams for one pion exchange between the two mesons bound in 
the molecule (A) and for the decay mode X — > D°D°ir (B). 



The effect of this reabsorption process can thus influence the partial width of the 
state in this decay mode. To obtain a qualitative estimate one can compare the average 
formation time Tf of the would-be-real-pion with the average time which elapses between 
the exchange of two subsequent pions T co u between the mesons. If Tf > T co u there is 
some probability that the pion is reabsorbed in the pion clouds which surrounds the D° 
meson, otherwise on average the reabsorption will not occur. This semiclassical analysis, 
albeit naive, is the only model independent possibility. 

The avreage formation time is proportional to the perpendicular component of the 
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pion four momentum k± in the rest frame of the X. From the two body decay kinematics 
we know that k± < fc*, the decay momentum in the D°* rest frame, thus: 

(68) T f ~ 1 

Since the D°* is not on shell fc* = ^X(q 2 , m 2 D , m 2 )/2y // q 2 . The off-shellness q 2 of the 
D°* meson is due to the exchange of pions with the D° inside the molecule. However, 
as the binding energy is extremely small the pions exchanged are soft, in the sense that 
a hard pion would break the molecule. 
In the soft-pion regime one finds that: 

(69) q 2 = {p - k) 2 = Mk_ m l + 2 ((k°f (fc 3 ) 2 - ki) 

where p is the four momentum of the two D-mesons. Furthermore we can neglect the 
temporal and longitudinal component of the pion four momentum: 

M 2 

(70) q 2 ~ -f- - m 2 D - 2ki => q 2 mm < q 2 < q 2 max 

where q 2 mm = M\j2 - m 2 D - 2(k™ ax ) 2 and q 2 nax = M 2 x /2 - m 2 D . When one takes 
into account the fact that Tf is a decreasing function of q 2 in the region of interest 
(q 2 > (mo + m^) 2 ) a range for the average formation time can be computed: 

(71) 1.5 x 1CT 23 sec < Tf < 1.8 x 1CT 23 sec 

which is of the same order of magnitude of the typical time of strong interactions. 
We now give an estimate of T co u in a scmiclassical approach. We describe the D°D°* 
5-wave molecule as two point-like particles which move in a one dimensional box of 
length ~ lOfm, the size of the molecule. The two particles exchange pions, i.e. interact, 
when their distance is smaller than 1 lfm, the range of strong interactions. To account 
for the quantum nature of the system we implement the uncertainty principle at each 
step, choosing the relative momentum from a Gaussian with mean value equal to the 
decay momentum fc* ~ 40 MeV and spread Afc ~ 1/Jf?, where R is the relative distance. 
The initial position of the particle is chosen from a Gaussian with zero mean value and 
AR = 10 fm. We put a cutoff on the relative momentum A — 200/300 MeV. The time 
evolution last for the average life-time of the X which is r = K/T(X) ~ 2.2 x 10 _22 sec. 
The evolution is repeated for several initial configurations. The outcome is stable against 
changes of the time step length and of the momentum cutoff and gives T co u ~ 10~ 22 sec. 
This means that T f < T coU and thus one expects that Y(X -> D°L»°7r ) ~ T(D°* -> 
D°ir a ). Any deviation from this prediction would seem unjustified in the molecular 
picture. 

We finally consider radiative decays. To describe the decays of the X in the tetraquark 
picture we need to introduce three amplitudes. The decay of a diquark-antidiquark bound 
state into a pair of mesons can occur through the exchange of a quark and an antiquark 
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belonging respectively to the diquark and the antidiquark. There arc indeed three differ- 
ent flavor configurations: the exchange of two light quarks A ([cq] [cq] — > [eg] [eg]) = A±, 
the exchange of two heavy quarks .4 ([eg] [eg] — > [eg] [eg]) = A 2 and the exchange of a 
light quark and a heavy quark .4 ([eg] [eg] — > [gg] [cc]) = A 3 . A\ and A2 account for 
X -> D°D *, while A 3 accounts for X -> J/ijmi:. 



X A 3 



X A 3 



(a) X J/t/ip (b) X -> J/t/iy 

Figure 22.: Radiative and hadronic decay of the X(3872) described with the same contact 
vertex A 3 . The radiative decay proceeds indeed through the hadronic transition X — > 
J/iPp. 



It has been experimentally confirmed that the pion pair in the decay into J/ifiir + TT 
comes mainly from a p. Thus we can write the partial width for this decay as: 



(72) 
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Since the amplitude for the radiative decay X — > J/ipj proceeds through the annihilation 
of a pair of light quarks into a photon, the hadronic part of the amplitude is the same 
as in the decay X — > J/ipTnr, Fig. 22. Exploiting the vector meson dominance one can 
write: 



(73) (J/ipj\X) = ( 7 |p) - J ( J /MX) = ^A 3 

7n p m z p 

The partial decay width can thus be written as: 

(74) hM . V 2Mx 
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Using f p = 0.152 GcV 2 [201] one obtains T(X J/tp^)/T(X -> J/^tttt) - 0.84, 
which is in rather a good agreement with the experimental value reported in Tab. XII: 
B(X -> J/ip^)/B{X -> J/tttt) ~ 0.3 ± 0.1. 

We can further exploit the result obtained for the width of X — > J/tpj in order to 
give an estimate of the decay width into J/i/> 77 . We can indeed compute the transition 
matrix element in terms of A 3 exploiting the coupling of the J/ip to the photon: 



(75) ( 77 |A) = ( 7 |J/V)^<JAH*) = M^{J/^i\X) = %^4^3 

The partial decay width can thus be written as: 



(76) r ( x^ 77) =wfe^ (AY i.vwsm 



3 H/J W7 8^M| 2M X 



Using = 1.254 GeV 2 [202] one obtains T(X -> 77 )/r(A -> J/^tttt) - 0.02, which 
has to be compared with the experimental value in Tab. XII: B(X — > 77 )/S(A — > 
J/tttt) < 0.01. The inconsistence of the theoretical prediction with respect to data is not 
dramatic if one takes into account the very strong assumptions made to derive Eq. (76). 

4'6.2. The 3940 family and other C = + states. A(3940) and X(4160) 
The possible assignment of V(3940) with the charmonium r] c (3S) with quantum numbers 
J PC = h would be supported by the following facts: 

• the non observation of the DD decay, 

• the lower states T] C (1S) and r] c (2S) are also produced through double charmonium 
production, 

• the predicted width for a state with a mass of 3943 MeV is - 50 MeV [203], 
in reasonable agreement sith the observed width. 

However the predicted mass is below the potential model estimate of <~ 4050 [204]. 
Moreover the observation of the A (4160), compatible with J PC = h and at the same 
distance from the predicted mass of 3 1 S as the A(3940), makes the picture more com- 
plicated. In fact the next radial charmonium excitation rj c (4:S) is predicted to be close to 
4400 MeV, so that it seems unlikely to accomodate both A (3940 and A (4160) as charmo- 
nium states. Additional important tests would be the study of the angular distribution 
of the DD* final state and the possible observation in 77 — > DD* [205]. 

Another interpretation which can not be excluded apriori is that of the A (3940) as 
a h hybrid charmonium, which can not decay to DD [74]. Spin dependent splittings 
place two states, h and 1 h , below the vector 1 with equal mass gaps of the order 
of 100 MeV [87, 88, 206, 207]. If the 1 — is identified with the y(4260), then it is not 
excluded that A (3940) is its lighter h partner. 
y(3940) and V(4140) 

The molecular nature of the V(3940) has been considered as a possibility by a number 
of authors [208, 147, 209, 55]. According to this picture, the hadronic wave function 
would be -^(\D* + D*~ > +\D*°D*° >), with J PC = 0++ or 2++. In this scenario 
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the Y(3940) would have a molecular partner, the Y(4140), with composition D* S D*, 
quantum numbers J PC = ++ or 2 ++ and a similar binding energy of about 80 MeV 
[210, 60, 211]. This idea is supported by the fact that the mass difference between these 
two mesons is approximately the same as the mass difference between the </> and u> mesons: 
TO Y(4i40)~ m Y(3940) ~ m^ — m^ ~ 210 MeV. However, with a meson exchange mechanism 
to bind the two charmed mesons, it seems natural to expect a more deeply bound system 
in the case that pions can be exchanged between the two charmed mesons, as in the 
D*D* , than when only 77 and <j> mesons can be exchanged, as in the D* S D* S system [12]. 
The molecule picture predicts that decays proceed via rescattering with decays to hidden 
and open charm states equally probable, so that decays to DD and DD* are foreseen. 
Another prediction of the molecular hypothesis is that the constituent mesons can decay 
independently, leading to D* s +Djj and D* s ~D+j [208, 39, 212]. Predictions for the 
radiative decays Y(3940)/V(4140) — > 77 in Ref. [55, 213] yield similar results for the 
jPC q++ an j 2++ assignments. In addition, a D* + D* s ~ molecule is predicted with mass 
of about 4040 MeV, decaying to J/ipp [61]. 

The conventional cc assignment of the Y(4140) is very unlikely since, being its mass 
above open charm threshold, it should have a large width. It was also shown in [214] 
that the Y(4140) probably can not be the second radial excitation of any of the P-wave 
charmonium states: x c j (J = 0, 1)- If it were the case, the branching ratio of the 
hidden charm decay, V(4140) — > J/i[)(f>, would be much smaller than the experimental 
observation [215]. 

The [cs][cs] tetraquark hypothesis for the V(4140) has been investigated by several 
authors [61, 216, 208]. However, in this scenario, it would decay with similar widths, of 
the order of 100 MeV, to hidden and open charm final states. 

With the present experimental information, also other explanations for the Y(4140) 
can not be excluded, such as hybrid charmonium [66, 61, 217] (its mass lyes in the range 
predicted by QCD for hybrids) or rescattering of D S D* S [218] or the opening up of a new 
final state channel [219]. Ref. [220] performed a study of the vector vector interaction in 
the framework of the hidden gauge formalism, finding three resonances with poles close 
to the masses of V(3940), Z(3930) and A" (4160) and quantum numbers J PC = 0++, 
2 ++ and 2 ++ respectively. In Ref. [221] their radiative open charm decays have been 
investigated. 
Af4350) 

The possible quantum numbers for a state decaying into J/ip(f> are J PC = 0++, 1-+ 
and 2 ++ , with the second possibility being exotic. Even with the modest information 
available for this state, different interpretations already appeared in the literature: a 
[cs][cs] tetraquark with J PC = 2++ [216], a molecular D*+D* s q state [222, 223], an 
excited P wave charmonium state S" 2 [224] and a mixed charmonium- D*D* state [225]. 
In Ref. [226] a QCD sum rules study has been performed to test whether the X(4350) can 
be an exotic J PC = 1 h D*D* (or D* Dq) molecular state. The mass value obtained is 
5.5 ± 0.19 GeV (4.92 ± 0.08), thus inconsistent with the experimental mass value. 
Z(3930) 

The measured properties of Z(3930) are consistent with expectations for the previously 
unseen l?Pi charmonium state, x' C 2 [203, 163], for which the predicted mass and width 
are 3972 MeV and 28.6 MeV (assuming the observed mass), respectively [72, 204, 203]. 
Also the measured two-photon production rate is consistent with the one of the x' C 2 [227]. 
Additional confirmations of this interpretation would come from the observation of the 
Z(3930) — > DD* decay with a branching ratio ~ 25% and the from the radiative decay 
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Z(3930) ^(2S)j with a partial width of the order of 100 KeV [204, 203, 205]. 

4'6.3. The 1~ family. The three 1~ states, V(4260), V(4360) and V(4660) are 
characterized by large total widths and by charmonium decay modes. Since their masses 
are higher than the D^D^ threshold, if they were 1 charmonium states they should 
decay mainly to D^D^*\ However, the observed Y states do not match the peaks in 
cross sections. Furthermore, the first available charmonium state is 
fy(3D), for which the quark models predict a mass of about 4500 MeV, which does not 
fit any of the new states. Therefore, the masses and widths of these three new Y states 
seem to be inconsistent with any of the 1 cc states [228, 229, 230]. An interesting 
interpretation is that the V(4260) is a charmonium hybrid [231, 232, 233]. Its mass lies 
infact in the ball park predicted by lattice QCD[234] and flux tube model. Actually, 
recent lattice simulations [235] and QCD string models calculations [236] predict that 
the lightest charmonium hybrid has a mass of about 4400 MeV, which is closer to the 
mass of the V(4360), thus not excluding the interpretation of the latter as an hybrid. 
A prediction of the hybrid hypothesis is that the dominant open charm decay mode 
would be a meson pair with one S'-wave D meson (D, D* , D s , D'*) and one P-wave D 
meson {D\, D s \) [237, 231, 99]. In the case of the V(4260) this suggests dominance of 
the decay mode DD\. Therefore, a large DD\ signal could be understood as a strong 
evidence in favor of the hybrid interpretation for the V(4260). Up to now the Belle 
experiment did not find evidence for such a signal [238, 239]. In the case of the F(4360) 
and V(4660), since their masses are well above the DD\ threshold, their decay rates 
into DDi should be very large if they were charmonium hybrids. Another prediction 
of the hybrid scenario is the existence of partner states: the flux tube model predicts a 
multiplct of states nearby in mass with conventional quantum numbers (0 h , 1 H , 2 h , 
1 ++ , 1 ) and states with exotic quantum numbers (0 H , 1 h , 2 H ) [240]. 

In Ref. [22] a possible interpretation of the y(4260) and V(4360) as D 1 D* molecular 
states has been proposed: two S'-wave mesons could be bound via pion exchange. A 
distinctive decay channel would be DD3ir. Further molecular interpretations of the 
y(4260) bound by meson exchange has been considered in Refs. [62, 241]. 

As already stressed, a crucial information for understanding the structure of these 
states is whether the pion pair comes from a resonance state. From the di-pion invariant 
mass spectra, there is some indication that only the V(4660) has a well defined interme- 
diate state consistent with / (980). Due to this fact and the proximity of the mass of the 
ip' - /o(980) system with the mass of the V(4660) state, in Ref. [82], the V(4660) was 
considered as a / (980) ip' bound state. If this interpretation of the V(4660) is correct, 
heavy quark spin symmetry implies that there should be a rf c — / (980) bound state 
[242]. This state would decay mainly into ri^im, and the authors of Ref. [242] predicted 
the mass of such a state to be 4616lg MeV. The enhancement at M — 4630 MeV in 
the A+A~ distribution has been suggested to be a FSI effect of the V(4660) — > A+A~ 
[243]. The V(4660) was also suggested to be a baryonium state [244], a canonical 5 3 Si 
cc state [245], and a tctraquark with a [cs] -scalar-diquark and a [cs]-scalar-antidiquark 
in a 2P-wave state [246]. 

In the case of V(4260), in Ref. [247] it was considered as a [cs]-scalar-diquark [sc]- 
scalar-antidiquark in a P-wave state. Studying the uncertainty in the determination of 
the orbital term, a mass of M = (4330 ± 70 MeV) was estimated, in nice agreement with 
the mass of V(4260) but also consistent with the mass of V(4360). However, from the rnr 
mass distribution, none of these two states, V(4260) and V(4360) has a decay with an 
intermediate state consistent with / (980) and, therefore, it is not clear that they should 
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have an ss pair in their structure. Besides, in Ref. [246] the authors estimate that the 
mass of a [cs]s = o[cs]s=o tetraquark in a P-wave state would be 200 MeV higher than the 
y(4260) mass. The authors of Ref. [246] found that a more natural interpretation for the 
y(4260) would be a [cg]s = o[«7]s=o tetraquark in a P-wave state. Other interpretations 
for the F(4260) appeared in the literature: a baryonium A c — A c state [248]; an S- 
wave threshold effect [71]; a resonance due to the interaction between the three, J/i/jtttt 
and J/ipKK, mesons [77]; a 4S 1 charmonium state [249]; an S- wave molecule p°Xa 
[250] or coxa [75]; an enhancement (not a true resonance) connected with the opening 
of the D*D* threshold and the coupling to the J/f/>/o(980) and J/-0cr(6OO) channels 
[251, 252, 253] The three Y states were also interpreted as non-resonant manifestations 
of the Regge zeros [254] . 

4'6.4. Charged States. As already stressed, the real turning point in the discussion 
about the structure of the new observed charmonium states was the observation by Belle 
Collaboration of a charged state decaying into iP'tt + , produced in B + — > K^'tt + [154]. 
Since the minimal quark content of this state is ccud, this state is a prime candidate for 
a multiquark meson. 

Z+(4430) 

There are many theoretical interpretations for the Z + (4430) structure. Since its mass 
is close to the D* D\ threshold, Rosner [66] suggested that it is an S'-wave threshold 
effect, while others considered it to be a strong candidate for a D*D± molecular state 
[64, 255, 67, 68, 69, 256]. 

Considering the Z + (4430) as a loosely bound S'-wave D*D\ molecular state, the 
allowed angular momentum and parity are J p = 0~ , 1 _ , 2 _ , although the 2~ assignment 
is probably suppressed in the B + — > Z + K decay by the small phase space. Among the 
remaining possible 0~ and 1~ states, the former will be more stable as the latter can 
also decay to DD\ in S-wave. Moreover, one expects a bigger mass for the J p = 1~ 
state as compared to a J p = 0~ state. The molecule explanation predicts that the 
Z + (4430) decays into D*D*n [64] through the decay of its consituents and tl){2S)iT via 
rescattering [240]. 

There is also a quenched lattice QCD calculation that finds attractive interaction 
for the D*Di system in the J p = 0~ channel [257]. The authors of Ref. [257] also 
find positive scattering lenght. Based on these findings, they conclude that althoug the 
interaction between the two charmed mesons is attractive in this channel, it is unlikely 
that they can form a genuine bound state right below the threshold. 

Other possible interpretations are tetraquark state [247, 258, 259, 260], a cusp in 
the D*D\ channel [261], a baryonium state [244], a radially excited cs state [262], or a 
hadro-charmonium state [263]. The tetraquark hypothesis implies that the £?+(4430) will 
have neutral partners decaying into ip'w /r) or r) c (2S)p° /u>. According to the tetraquark 
picture the decay channels should be DD* , D*D* , J/ipir, J/ipp, VcP and ip(2S)w, but not 
DD because of its spin-parity [65] Further discussions about its production and decay 
can be found in Refs. [66, 264, 265, 266] and [267], respectively. The tetraquark model 
also predicts a second nearby state with mass ~ 4340 MeV, decaying into i/>'7r + [268] . 

£+(4050) and Z^(4250) 



Due to the closeness of the Z+(4050) and Z+(4250) masses to the P*P*(4020) and 
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DiZ)(4285) thresholds, these states could be interpreted as molecular states or threshold 
effects. However, since the mass of ^(4050) is above D* D* threshold, the molecular 
interpretation is disfavoured, even if studies present in the literature give contradictory 
results: in Ref. [76], using a meson exchange model, strong attraction for the D* D* 
system with J p = + is found, while using a boson exchange model, the author of 
Ref. [56] concluded that the interpretation of Z^~(4050) as a D*D* molecule is not favored. 
In the case of ^(4250), using a meson exchange model, it was shown in Ref. [62] that 
its interpretation as a D\D or D^D* molecule is disfavored. 

4'7. Summary and outlook. - A large number of new Charmonium states have been 
measured and the knowledge of their characteristics has significantly improved. Fig. 23 
shows the observed states overlaid on the regular charmonium states as predicted by 
the potential models [2]. The most likely J CP assignment is shown. Charged states are 
labelled with the most likely J p assignment. Colours identify the grouping of the states 
within the previous paragraphs: red for the individual states (Sees. 4T and 4'3), green 
for the 3940 family (Sees. 4"2), blue for the states around M = 4140 MeV (Sec. 4'3), 
purple for the 1 statesSec. 4'4), and orange for the charged states (Sec. 4'5), Black 
states are traditionally considered as regular charmonium, although ^(4040) seems to 
behave as an exotic particle since it does not match any potential model and is close in 
mass to a charged state with the same J p . In summary of the observed states above open 
charm threshold only the X(4160), Z(3940), F(3940), and V(3770) are good candidates 
for regular charmonium: n c (3S), Xc,2(2-P), Xc,o(2-P), and ip(lD) respectively. 
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Figure 23.: Observed states with hidden charm above the open charm threshold. The 
theory predictions are according to the potential models described in Ref. [2]. 



Besides observing an increasing number of exotic states we are still far from hav- 
ing clear assignments for each of the states between different possible interpretations. 
Moreover, a new spectroscopy implies the existence of a large number of states, whose 
absence would have to be justified. A perception of the status of the global picture and 
an indication on where to search can be obtained from a comparison between possible 
spectra and observed states. 

Expected tetraquark spectra were derived, under assumptions which lead to uncer- 
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tainties 0(100 MeV), in Sec. 2'1.1. Fig. 24 and Fig. 25 show a comparison between 
the expected and observed spectra for the first two radial excitations of J = + ,l~, 
and 1 + spectra and the observed states. Besides the V(3872) which is assumed to be 
a tetraquark when building the model and that would naturally be constituted by two 
states close in mass, the following states have a match within the 100 MeV : V(4350), 
Z(4430), X(4160), y(4260), V(4350) and V(4660). Such a small number of matched 
states is opposed to the large number of needed states: 18 (27) for each of the J = 0(1) 
multiplets. While it can be argued that production and decay mechanism can differen- 
tiate between the states and that experimental sensitivity differs significantly between 
final states (for instance D s mesons are much more difficult to detect than D mesons), 
one striking observation is that no attempt is made to search for the strange states, that 
would decay into charmonium plus a kaon. 

Predictions in the case of molecules are more difficult. It is easy to classify the 
masses around which a molecule could possibly lay by computing all sums of the masses 
of two mesons with correct quantum numbers (see Fig. 26). Here all pairs of either 
a charmonium and a light meson, neglecting the scalar nonet which is here treated as 
a tetraquark itself and states with a width larger than 50 MeV, or two open charm 
mesons. Each case should be considered separately to estimate the production cross 
sections and the binding energies. Since the molecules have masses lower than the sum 
of the constituent mesons it is interesting to search for their decays in the final states at 
lower masses. 

Among the observed state the V(3872) is the only state that matches a threshold 
within 10 MeV. 

5. Bottomonium 

The discovery of exotic charmonium resonances that do not fit in the standard charmo- 
nium scheme suggests the presence of their bottom companions, with similar properties. 

The interest on this field has been boosted by the Belle discovery of an anomalously 
large Y(nS)n + n~ (n = 1, 2) production around the T(5S) resonance [269]. Since the 
T(nS)Tr + Tr~ channel is the equivalent, in the bottom sector, of the ip(nS)ir + Tr~ channels 
preferred by some of the exotic charmonia, the interpretation of this result as the effect 
of an exotic bottomonium with a mass around the T(55) peak is quite natural ( x ). 

The current predictions for tetraquarks with bottom content have been already pre- 
sented in Sec. 2T. Here we report the searches performed by Belle and BaBar (Sec. 5'1) 
and review some interpretations of the experimental results (Sec. 5'2). 

5T. The search for exotic bottomonium. - Searches for exotic bottomonia have been 
performed at the i?-Factories exploiting inclusive and exclusive techniques. In the in- 
clusive searches, structures are searched in the e + e~ — > hadrons cross section, looking 
in particular at the hadronic ratio R = er(e + e~ — » hadrons) / a (e + e~ — > ,u + /i~). In the 
exclusive analysis, specific final state channels are searched for, with a particular atten- 
tion devoted to the bottom equivalents of the decay channels preferred by the exotic 
charmonia, in particular T (nS)ir + n~ . 



( 1 ) One could also hypothesize the presence of an exotic component mixed with the T(5S) to 
form a single resonance. 
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Figure 24.: Expected tetraquark spectra as derived in Sec. 2' 1.1 for the first and second 
radial excitations of states with J p = + , 1 + with the observed states superimposed. Color 
coding is the following: [cl][cl], where I = u,d are purple triangles (charged), light green 
double crosses (neutral with C = +), and light blue crosses (neutral with C = —); [cl][cs\ 
and charge conjugates are red boxes; [as] [cs] are dark green double squares(C = +), and 
dark blue squares (C = —). 



5*1.1. Spectroscopy from the Belle T(55) campaign. We already mentioned the obser- 
vation, made by the Belle collaboration, of an anomalously large T(nS')7r + 7r _ around the 
T(55) [269]. This analysis is based on the first data (21.7 fb _1 ) collected at the center of 
mass energy of the T(55) resonance. The authors look for the T(rtS')7r + 7r~ final states, 
possibly produced in association with an ISR photon, by selecting events with four tracks, 
two of which are identified as muons and have an invariant mass consistent with the 
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Figure 25.: Expected tetraquark spectra as derived in Sec. 2 '1.1 for the first and second 
radial excitations of states with J p = l - with the observed states superimposed. Color 
coding is the following: [cl][cl], where I = u,d are blue crosses, [cl][cs\ and charge conju- 
gates red boxes, and [cs] [cs] are purple squares. Different values of C are degenerate in 
the model and therefore they are not separated in the plot. 



T(nS) — > fi~~ hypothesis. If a resonance of mass M decays into T(nS)7r + 7r _ , the dis- 
tribution of the invariant mass difference AM = M (ix + pT n + tt~ ) — M(n + fi~) is expected 
to peak at AM = M — Mx(nS)- A sharp peak corresponding to M ~ Mx<5S) ~ 10.865 
GeV is found in both T(lS')7r + 7r _ and T(25)7r + 7r _ channels, which is inconsistent with 
the naive expectations from the scaling of the T(4S I ) rates [270], as shown in Fig. 27. 

Subsequently, the Belle collaboration performed an energy scan with y/s between 
10.83 and 11.02 GeV, in six steps, for a total integrated luminosity of ~ 7 fb^ 1 [271]. 
In this case, the authors look for events around AM — y/s — My/ n s) (i.e. without ISR 
production), in order to measure the T(2S')7r + 7r _ rate in each of the six energy points 
and determine the resonance shape. The analysis confirmed the large rate observed in the 
previous test, adding also the observation of the T(3S')7r + 7r~ channel. The combined fit of 
the resonance shapes with a single Breit-Wigner PDF, interfering with a flat continuum, 
gives a mass of 10.8884±°;°°^±o.0012 GeV and a width of 30.7±^±3.1 MeV, that differ 
by about 3cr and 5a respectively from the values quoted for the T(55*) in the PDG [272]. 
This observation triggerd the hypothesis of two near but separated structures around 
10.88 GeV. 

5T.2. The BaBar scan above the T(4S I ). Before the works by Ali et al., Ebert et al. 
and Wang et al. were published, a naive scaling of the exotic charmonium masses to 
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Figure 26.: Sum of the masses of any pair made by either a charmonium and a light 
meson (blue dotted lines) or two open charm mesons (orange full line) for C = +(top) 
and C = — (bottom). The assigned J PC are computed assuming L = and null charge. 
Black dashed lines represent the observed exotic charmonium states. 



the bottom case suggested the presence of exotic bottomonia in the region between the 
T(45) (10.58 GeV) and the T(65) candidate (~ 11 GeV). In particular, it is well known 
that the masses for the regular bottomonia can be obtained, at a good approximation, 
by scaling the masses of the regular charmonia by the mass difference between the T(1S) 
and the J I if). If it is assumed to work also for the exotic mesons, the masses of the 
y(4260), F(4350) and F(4660) can be scaled to "predict" three exotic 1 bottomonia 
at 10.62, 10.71 and 11.02 GeV, that could be directly produced in e + e~ collisions. The 
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Figure 27.: Distribution of AM = M (/i + [i~ ir + ir~ ) — M(fi + /.i~) obteined at the center 
of mass energy of the T(5S I ) (top), compared with the distribution of the same variable 
at the T(4S*) (bottom), where the three arrows indicate, from the left, the T(25), T(3S) 
and T(45) peak. 



BaBar collaboration decided to investigate this possibility by performing an energy scan 
between 10.54 and 11.20 GeV. Steps of 5 MeV were performed, by collecting about 
25 pb^ 1 per step, for a total of 3.3 fb~ . An additional data set of 600 fb _1 was also 
collected in 8 steps around the T(6S) candidate (10.96 to 11.10 GeV). This scan overtakes 
by a factor 30 in statistics and a factor 4 in step fineness the previous one, performed at 
Cornell [273, 274]. The search for exotic bottomonia can be performed in these data set 
exploiting both inclusive and exclusive reconstruction. 

At present, only the results of the inclusive analysis have been published [275]. 
Hadronic final states are selected by requiring at least three tracks, with a reconstructed 
energy greater than 4.5 GeV, while event shape variables are used to reject the QED back- 
ground. Di-muon events are selected by looking or two tracks, collinear within 10°, with 
an invariant mass of at least 7.5 GeV. The bb (wy) contribution Ri, to the total hadronic 
ratio is measured, by subtracting from the total hadronic yield the contribution from 
continuum (e + e~ — > qq, q — u,d,s,c) and two-photon (e + e~ — > e + e _ 7*7* — > e + e~Xh) 
events, estimated below the BB threshold and scaled according to the expected y/s de- 
pendence, while the ISR (e+e~ — s- Y(nS)~fisii) contribution is included in Rt. The 
result of the analysis is showed in Fig. 28. The authors also performed a simple fit in the 
region between 10.80 and 11.20 GeV, including two Breit-Wigner resonances and a flat 
continuum background, allowing a partial interference between the three components: 
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a = \A nr | 2 + \A r + A W860 BW(M W86 o, Tioseo) + ^1102051^(^^11020, rno2o)| 2 - The mea- 
sured T(55) mass and width are 10.876 ± 0.002 GeV and 43 ± 4 MeV respectively, but 
are affected by large systematic uncertainties related to choice of the fit model. Recently, 
the Belle collaboration updated the analysis of [271] , by including a measurement of Rb 
and a fit that replicates the BaBar fit and confirms these results. At present, the Belle's 
analyses provide a discrepancy of only 2.2cr (corresponding to 9 MeV) and 1.4ct in mass 
and width between the inclusive and the exclusive measurement. 

5'2. Interpretations of bottomonium results. - The interpretation of the results ob- 
tained so far in the search for exotic bottomonia is still quite controversial. On one side, 
the observation of a large Y(nS)ir + ir~ rate around the T(55) is in contrast with a naive 
rescaling of the corresponding T(4S*) rate, and can suggest the presence of an exotic 
resonance lying near the T(55); this hypothesis is also strengthened by the study of the 
invariant dipion mass spectra and the angular distributions of the observed decay [276] . 
On the other hand, although a well established model for this decay of the standard bot- 
tomonium is still missing, some works have shown that it is possible in principle to bring 
back the Belle observation into the standard scheme [277, 278]. In the Meng and Chao's 
work a rescattering model is used, where a T(55) first decays into a BB pair, then a B 
meson exchange produces a lower T and a scalar resonance (a or / ) decaying into n + n~ . 
Instead, in the Simonov and Veselov's paper, an effective Lagrangian for qq, qqir and qqinr 
production is introduced. In both cases, parameters are fixed independently of the T(55) 
measurements (from charmonium or lower T's) and a large T(55) — )■ Y(nS)ir + ir~ rate 
is predicted, in rough agreement with the experimental results. 

The observation of a significant difference between the measured mass and width of 
the T(5S) in inclusive and exclusive decays is also not surprising even in the standard 
framework. A proper description of the hadronic ratio R near the open charm and bottom 
thresholds requires to take into account the mixing of the pure T(nS) states, introduced 
by their coupling with the OZI-allowed states £>(*)£>(*) and £(*)#(*) [279, 280]. Such a 
kind of treatment, for instance, allowed to predict the two structures observed by BaBar 
between 10.60 and 10.75 GeV without introducing any exotic resonance. An important 
conclusion of these studies is that the total hadronic cross section is not expected to be 
a simple superposition of Breit-Wigner resonances, and significant differences can arise 
between the inclusive and exclusive shapes of the cross section. Hence simple fits based 
on a superposition of Breit-Wigner resonances are not well grounded and the results 
are difficult to interpret, and in absence of a consistent treatment of the inclusive and 
exclusive processes within the same theoretical framework, no conclusion can be drawn 
from the small differences observed in the shapes. 

Recent efforts to better understand the hadronic ratio, in particular in the region 
between the T(55) and the T(65), have been also performed by Ali et al. [281] and by 
van Bcvcren and Rupp [282]. In particular, in the Ali et al. paper the exotic resonance 
responsible for the large T(n«S')7r + 7r _ production is identified with a possible peak noticed 
in the BaBar Rb measurement around 10.91 GeV, although from an experimental point 
of view it doesn't seem to have enough statistical significance to claim an observation. 
Instead, in van Beveren and Rupp's paper, a quite uncommon interpretation of the BaBar 
result is presented: the peak at 10.58 GeV, usually identified with the T(45) resonance, 
is interpreted as a threshold effect, while the T(45) is found at 10.735 GeV; moreover, 
the peaks around 10.88 and 11.02 GeV, usually assigned to the T(55) and T(6S), are 
identified with the T(3D) and the T(55), respectively. 

Finally, in Fig. 28 we illustrate the BB thresholds in whose neighborhood molecular 
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Figure 28.: BB thresholds (left) and expected tetraquark 1 states (right), superim- 
posed to the Rb spectrum measured by the BaBar collaboration. 



states could be expected, and the tetraquark states predicted in this work (see Fig. 5), 
superimposed to the Rb spectrum measured by the BaBar collaboration. 

6. Outlook and conclusions 

It is now almost seven years since the first heavy quarkonium exotic states has been 
observed and the number of exotic states has significantly increased. In order to explain 
all these states one needs to go beyond the assumption that mesons can only be made of 
a quark and an anti-quark and consider the possibility to be observing bound states of 
four quarks or two quarks and valence gluons. These novel states of aggregation are not 
a new idea and it was indeed already in the foundations of the quark model [5]. Also, 
the hypothesis that the mesons of the scalar nonet are actually tetraquarks can be dated 
back to the '70s [283]. But only now, with the presence of several exotic candidates with 
heavy quarks the possibily to have a global picture is becoming real. 

Unfortunately, such a spectroscopy made of a larger number of constituents (four 
quarks for tetraquarks and molecules and two quarks and gluons for hybrids) is much 
reacher than standard spectroscopy and therefore a significantly larger number of states 
could potentially exist. To the aim of building this global picture, a complete set of 
theoretical predictions and a systematic experimental search is needed. The work in 
the theoretical path is progressing although strong interactions are extremely hard to 
compute, even in presence of heavy quarks, and therefore predictions are extremely hard 
and uncertain. From the experimental point of view the amount of work ahead is if 
possible even larger. 

The status of the experimental observations is extremely fragmented. Pictorial sum- 
maries, separated by production mechanisms, are in Figs. 29-31. The exotic states have 
been observed always in only one production mechanism and often in only one final state. 
None of the states has been searched systematically in all final states. Sometimes the 
analysis of a given final state for a given production mechanism is missing, but often it 
has either been performed only in a limited mass range or the invariant mass spectrum 
has been published without a fit to the possible new state. This is mostly due to the 
fact that some of the analyses did not show a significant signal themselves and they were 
published before a new state was observed. As an example, the invariant mass spectra of 
a charmonium and a photon in Fig. 9 are focussed on the A(3872) region and no infor- 
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mation is available outside it. On the other side in the same Fig. 9, the J/iprj invariant 
mass spectrum is published, but not fitted for all the possible new states. Whatever 
the cause is, the critical point is that building a global picture requires having either an 
observation or a limit in each final state for each candidate new state: limits allow to 
quantify the level at which a decay mode has not been observed and can show whether a 
signal could have been observed in a final state or the efficiency, the branching fraction 
or the background level would make it unobservable. It is also worth remarking that 
although finding new decay modes of these states might not necessarily cast light on the 
nature of these particles, they would in any case concur to the evidence of the existence 
of the states that suffer from lack of statistics. 

Looking into this "observational" tables in detail, B decays (top of Fig. 29) are the 
most studied, but there is a significant amount of missing fits ("M/F"). Particularly 
severe is the lack of analysis of the baryonic spectra {pp, shown in Fig. 19, and AA) 
since baryonic decays are a signature of tetraquark states. Some other modes have never 
been studied, mostly because the number of expected events is very low. Nonetheless 
the study of B — > ip(2S)nirK decays should be relatively clean, while D*D* and above 
all D^Di*^ suffer from the low branching fractions of the observed states. 

States produced in conjunction with an initial state radiation (ISR) photon (bottom 
of Fig. 29) have an unambigous J CP assignment and arc therefore less analyses are 
needed to establish their properties. Nonetheless it is striking to see that a large fraction 
of analyses have been carried out exclusively for the first observed exotic state, the 
y(4260). It can also be noticed that no search is published involving Di*^ mesons: while 
the efficiency is expected to be very low, background should be low as well and surprises 
can always arise. Finally the F(4660) has been object of one of the combined analyses we 
are advocating here [153]: two states apparently different, observed in *l)(2S) and A C A C 
final states, if fitted under the same ansatz were found to be consistent with being the 
same and interesting ratios of branching fractions were measured. 

On the recoil of a J/tp and in 77 interactions ( Fig. 30 ) only C — + neutral states 
can be observed. This restricts the number of final states of interest. Also, the low 
multiplicity of these decays and the large missing momentum in the case of 77 decays 
makes these analyses experimentally challenging. On the other side C = + states are the 
least known and reinforcing the evidence of the signals would help. It is also interesting 
to notice that, mostly due to statistics, the recoil to any other particle but the J/ip has 
not been investigated. Given the selection rules the recoil to Xco or Xc2 would be very 
interesting. 

Searches at hadron colliders (top of Fig. 31) have the advantage to have very large 
samples. On the other side backgrounds are very high and therefore final states with 
too high multiplicity and above all with neutral particles are not at reach. Finally, 
it is hard to extract information on production cross sections and therefore branching 
fractions. Nonetheless, it is clear that a systematic search in pp collisions is likely to 
clarify the picture significantly and it is a pity that it is still missing. This holds also and 
in particular for the charged states - for instance the longly debated Z + (4430) should 
give a large signal in the J/ipir + spectrum - and in the bottomonium - where no other 
experiment is in the position of studying inclusive Y(nS)ir + Tr~ and T(nS')7r ± spectra. 

Concerning the searches of charged exotic states, the most striking signature of states 
made of more than two quarks, very few searches have been conducted in B decays. We 
believe that for each exotic neutral spectrum the corresponding charged state should be 
searched for completeness and, as shown in the bottom of Fig. 31 only five combinations 
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of final states and exotic states has been searched for. As an example, no information 
has been extracted from Fig. 18 on the charged partner of the A"(3872), Z(3870) in our 
tabic, which has long been pointed out as a critical state to search for. Moreover among 
the four quark bound states there must be states that contain a single s quark (see for 
instance Fig. 24 and Fig. 25), their mostly distinctive signature being a charmonium plus 
a charged kaon. These states could be searched in B decays in association with a ss state 
or inclusively in pp collisions. 

Concerning the bottomonium states, the searches are extremely limited and there no 
confirmed evidence of exotic states yet. The search potential of B-Factories is limited by 
the fact that bottomonium at masses higher than the open bottom threshold can only be 
produced in conjunction with initial state radiation and that in addition the accessible 
mass range is full of threshold openings (see Sec. 5'2). This restricts the search to 
J PC = 1 states with significant signatures in exclusive final states. Besides the lack of 
information about most of the possible final states, the real contribution could come from 
hadron colliders, where the search could be inclusive over Y(nS)X final states (where X 
is made of charged tracks) and extended to exotic states with open bottom, eventually 
even with open charm. 

In conclusion, seven years of discoveries of exotic particles need to be followed by a 
systematic study of the possible new spectroscopy to be able to give a global and definite 
picture. To this aim still a lot can be extracted from existing data of _B-Factories and 
Tevatron. Nonetheless in the present generation of experiments statistics is extremely 
low and we will not be able to have a convincing picture until the results of ultra-high 
intensity machines, LHC and the Super Flavour Factories (SuperB and/or SuperBelle) , 
will be available. 
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Table X.: Design parameters of the three B-Factories. 





CESR 


KEK-B 


PEP-II 






LER HER 


LER HER 


Energy (GeV) 


5.29 


3.5 8.0 


3.1 9.0 


Collision mode 


2 mrad 


11 mrad 


Head-on 


Circumference (m) 


768 


3018 


2199 


/3*/ /3* (cm) 


100/1.8 


100/1 100/1 


37.5/1.5 75/3 


£x/ £y 


0.03/0.06 


0.05/0.05 


0.03/0.03 


e t/ e y (7rrad — nm) 


210/1 


19/0.19 19/0.19 


64/2.6 48.2/1.9 


relative energy spread (1CP 4 ) 


6.0 


7.7 7.2 


9.5 6.1 


Total Current (A) 


0.34 


2.6 1.1 


2.14 0.98 


number of bunches 


45 


5120 


1658 


RF Frequency (MHz)/ Voltage (MV) 


500/5 


508/22 508/48 


476/9.5 476/17.5 


number of cavities 


4 


28 60 


10 20 



Table XL: Design parameters of BEPC, CESR-c and the upgraded version of BEPC. 





BEPC 


CESR-c 


BEPC-II 


Max. energy (GeV) 


2.2 


2.08 


2.3 


Collision mode 


Head-on 


±3.3 mrad 


22 mrad 


Circumference (m) 


240 


768 


240 


/3*/ 13* (cm) 


120/5 


94/1.2 


100/1.5 


&/Cy (10- 4 ) 


350/350 


420/280 


400/400 


ej/ej(7r rad — nm) 


660/28 


120/3.5 


144/2.2 


relative energy spread (10~ 4 ) 


5.8 


8.2 


5.2 


Total Current (A) 


0.04 


0.072 


0.91 


number of bunches 


1 


24 


93 


RF Frequency (MHz)/ Voltage (MV) 


200/0.6-1.6 


500/5 


500/1.5 


number of cavities 


4 


4 


2 
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Table XII.: Measured X(3872) product branching fractions (PBF), separated by pro- 
duction and decay mechanism. When more than a publication is present the combination 
is performed assuming gaussian uncorrelated errors. The last two columns report the 
results in terms of absolute X branching fraction (B/it) and of the branching fraction 
normalized to J/ipTrn (Rfu) as obtained from the global likelihood fit described in the 
text. Ranges and limits are provided at 68% and 90% C.L. respectively. 



B Decay mode 


X decay mode 


PBF(xl0 5 ) 


Bfu 


Rfu 


XK± 


X ->■ J/tpTVTY 


0.82±0.09 [130, 131] 


[0.035,0.075] 


N/A 


XK° 


X — > J/ipnn 


0.53±0.13 [130, 131] 


N/A 


N/A 


XK* 


X -+ D*°D° 


13±3 [127, 132] 


[0.54,0.8] 


[7.2, 16.2] 


XK° 


X D*°D° 


19±6 [127, 132] 


N/A 


N/A 


XK 


X -S- J/i/>7 


0.22±0.05 [133, 134] 


[0.0075,0.0195] 


[0.19,0.33] 


XK 


X -s- ip(2S)j 


1.0 ±0.3 [133] 


[0.03,0.09] 


[0.75,1.55] 


XK 


X — > 77 


<0.024 [135] 


< 0.0004 


< 0.0078 


XK 


x j/i>n 


< 0.77 [136] 


< 0.098 


< 1.9 


XK 


X -S- J/tPtvtytt" 


" [134] 


[0.015,0.075] 


[0.42,1.38] 


XK* 


X —¥ J/TplYTT 


< 0.34 [131] 


N/A 


N/A 



Table XIII. : Measured J PC ', masses, and widths of the "3940 family" of states. 



State 




Mass (MeV) 


Width (MeV) 


X(3940) [141] 
Y(3940) [Belle] [53] 
Y(3940)[BaBar] [54] 
Y(3915) [142] 
Z(3940) [143] 


±+ 
[0,1,2] ±+ 
[0,1,2] ±+ 
[0,1,2] ±+ 

2++ 


3942±l(stat.) ± 6(sys.) 
3943 ± 13 
39U.6±l l(stat.) ± 1.9(sys.) 

3915 ±3(stat.) ± 2(sys.) 
3926±2.7(stat.)±l.l(sys.) 


37+ll(stat.)±8(sys.) 

87 ± 22 
33t\ 2 (stat.) ±5(sys.) 
17±W(stat.)±3(sys.) 
21.3 ±6.8 ±3.6 
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